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Abstract 

We analyze the high energy scattering of hadrons in QCD in an 
effective theory model inspired from a gravity dual description. The 
nucleons are skyrmion-like solutions of a DBI action, and boosted 
nucleons give pions held shockwaves necessary for the saturation of 
the Froissart bound. Nuclei are analogs of Blon crystals, with the 
DBI skyrmions forming a fluid with a fixed inter-nucleon distance. 
In shockwave collisions one creates scalar (pion held) “fireballs” 
with horizons of nonzero temperature, whose scaling with mass 
we calculated. They are analogous to the hydrodynamic “dumb 
holes,” and their thermal horizons are places where the pion held 
becomes apparently singular. The information paradox becomes 
then a purely held theoretic phenomenon, not directly related to 
quantum gravity (except via AdS-CFT). 
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1 Introduction 


In a series of papers m ei ei la ei ei it was shown that the high energy, small angle (large 
s, fixed t) scattering in QCD can be described through a simple cut-off AdS gravity dual 
a la Polchinski-Strassler [7], where the AdS space ends at an IR brane. Above the gauge 
theory Planck scale Mp = Nc/ 4 Aqcd, the scattering is dominated by black hole creation in 
the gravity dual. The maximal Froissart regime in QCD 

a = — -ln 2 s (1.1) 

m z 

corresponds in the gravity dual to black hole creation that occurs effectively on the IR brane. 
Then in the gravity dual we have an effectively 4d scattering in a gravity theory with a mass 
(the KK mass obtained by reducing gravity onto the IR brane), reducing to gravitational 
shockwave scattering in the large s limit, and creating black holes. It was shown in |2j that 
this model exactly matches the 1952 Heisenberg model for the saturation of the Froissart 
bound 0, a model in which one analyzes pion held shockwave scattering. Moreover, one 
needs a nonlinear DBI-type action for the pion held to obtain saturation of the bound. 
Heisenberg takes the action 


S = T(3) J d4x \l 1 + + m2 ^ 2 ( L2 ) 

Because of the matching of the two descriptions of the Froissart saturation, in |T it was 
proposed that there should be a pion held “soliton” created in the QCD collisions, radiating 
thermally at a given temperature, and this object was identified with the hreball observed 
at RHIC (see also urn for a different view of how the RHIC fireballs are related to black 
holes, in J\f = 4 SYM). 

A different use of the pion held was put forward in the program of Skyrme-like models 
muni esi mi- There one tries to hnd the nucleon as a topological solution of an effective 
pion held action with higher derivatives. The nonlinear (or linear) sigma model action 
does not admit solitons, but if one adds a higher order correction solitons can be obtained 
mi in ns. The original Skyrme model had a particular correction, but one obtains the 
same results with a large class of corrections. In particular, Pavlovskii El has analyzed a 
DBI-like action, that reduces to the sigma model at low energies and contains an infinite 
series of higher order corrections, summarized in the DBI square root. The action 

s = \j d 4 xflA\ CD [\j 1 + ( e -ir7h^ e ^ T 7h) 2 /(A QCD ) 2 - 1] (1.3) 

admits Skyrme-likc topological solitons that could be identified with the nucleons. 

In this paper we will try to put these two approaches together, and hnd an effective held 
theory description that will encompass both the hadrons at rest -nucleons in a Skyrme-like 
model- and the high energy scattering of hadrons. For the gravity dual description of the 
latter we often found direct effective held theory interpretation, so now we would like to use 
the intuition gained from the gravity dual to set up a purely effective held theory description 
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of the scattering. We know that QCD itself should give a good description of the scattering, 
but a first principle QCD description is very hard, and the gravity dual description is suited 
to an effective held theory description anyway, so that’s what we will find. However, this 
will not be a usual effective held theory, since if the Skyrme-like soliton only needs a few 
correction terms to the action, the large s scattering needs an action valid at energies well 
above the natural cut-off, A qcd- 

We will find an action that admits both a Skyrme-like topological soliton of the type 
of Pavlovskii and solutions with horizons that radiate thermally at a given temperature, 
corresponding to the “fireballs” observed at RHIC. The hreball solutions will be very similar 
to the “dumb holes” of Unruh EE|> obtained in sonic booms. The hydrodynamic equations 
of the sonic boom are exactly analogous to the scalar equations of motion in our theory. 

The paper is organized as follows. In section 2 we will review the high energy (large s) 
scattering in gravity duals of QCD, in section 3 we will look at various DBI actions and their 
Blon solutions. In section 4 we will boost the Blon solutions and compare to the boosting 
in the gravity dual model for QCD scattering. In section 5 we will analyze fireball-like 
Blon solutions with horizons and compare them with the “dumb holes”, thus calculating 
their temperature. In section 6 we will look at SU(2) actions and topological solutions, 
in particular the Pavlovskii solution. Section 7 contains the bottom line of the paper. We 
present our proposed DBI action, argue for its form from the gravity dual perspective together 
with QCD arguments and find its “SkyrBIon” solution. In section 8 we look at high energy 
scattering in our model and in section 9 we conclude. 

2 High energy scattering in QCD gravity duals 

In this section we review the description of large s, fixed t scattering via gravity duals of 
QCD. 

Following Polchinski and Strassler [7J, the amplitude for scattering in QCD can be found 
from a gravity dual by multiplying with wave functions in the extra dimensions and inte¬ 
grating over these extra dimensions 

A(p) = / drd 5 Qy/gAstring(p) (2.1) 

^ i 

We do not know of course the gravity dual of QCD, but we know that it should look like 
AdS$ x X 5 space modified in the IR (and maybe in the UV). 

ds 2 = -=^dx 2 + —jrdv 2 + R 2 ds\ = e~ 2y/,R dx 2 + dy 2 + R 2 ds 2 x (2.2) 

rz r z 

Polchinski and Strassler proposed that one could obtain a lot of information just by putting 
a sharp cut-off in the IR (at r min ~ R 2 A). In the large s, fixed t regime, the scattering in 
the gravity dual is concentrated close to, but not on the IR cut-off (IR brane) j2J. 

But high s, fixed t scattering in a gravitational theory was shown by’t Hooft El to be 
described by scattering Aichelburg-Sexl gravitational shockwaves [TH 

ds 2 = 2 dx + dx~ + (dx + ) 2 ^(x i )6(x + ) + dx 2 (2.3) 
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where the function $ satisfies the Poisson equation 


A d _ 2 $(i ! ) = — 16ttGp5 d 2 (x l ) 


(2.4) 


He showed that for s < Mpi one can treat one particle as a shockwave and the other as a 
null probe (geodesic), and proposed that above M P i both particles should be shockwaves, 
and one will create black holes. This case of s Mpi was analyzed (and the cross section 
for black hole production was computed) in flat 4 dimensions in jl9| and extended to higher 
dimensions and curved space in I]. As suggested in on on general arguments, it was found 
that the cross section in flat space grows like a power law 

a ~ r 2 H ~ (2.5) 


One can find solutions for gravitational shockwaves in curved spaces of gravity dual type 
1U> an d we find that the solutions still look like (|2.3|) in the given background, where now $ 
satisfies the Poisson equation in the background. 

In [21 Ej it was found that the behaviour for the cross section for black hole formation in 
the gravity dual translates into the same kind of behaviour for the QCD cross section, and 
we have the following regimes. 

Above the Planck scale, which in gauge theories corresponds to Mp = Nr/ 4 Aqcd with 
Aqcd the scale of the lightest glueball excitation, and in real QCD would be about 1-2 GeV, 
in the gravity dual we form small black holes, that “feel” only flat space. Correspondingly, 
one finds that the shockwave profile $ is 


167 tGd 1 

Qd -3 {D — 4)r D-4 r D ~ 4 


( 2 . 6 ) 


By applying the formalism for shockwave scattering with black hole formation and translating 
to QCD with (12.111 one finds that 

a ~ (2.7) 

both in the gravity dual and in QCD. 

As one increases s beyond Ep = N 2 /R in the gravity dual and above Ep = N 2 Aqcd in 
the gauge theory, about 10 GeV in real QCD, the black holes produced in the gravity dual 
start to “feel” the curvature of space. One finds that for consistency, in the AdS^+i x Xj 
gravity dual Xj needs to be large (with scale much larger in the IR than that of AdS), and 
the shape of the shockwave profile becomes 


, K] Rjr i i 

(J) - rv./ -- - 

— l)+d /y*W 

and by applying the formalism one finds 


( 2 . 8 ) 


<7 ~ 


s 1/n = s 1 ' 11 


(2.9) 


both in the gravity dual and in QCD. 
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Finally, the last regime corresponds to the maximal Froissart behaviour. Above an un¬ 
known energy scale E F , that depends on the details of the gravity dual in the IR, but in real 
QCD the experimental data suggests it should be between 100 GeV and 1 TeV, the black 
holes created in the gravity dual are so large that they reach the IR cut-off and get stuck 
there (since the scattering in the dual happens mostly near the IR). Thus the gravity dual 
scattering effectively happens on the 4d IR brane, and creates 4d black holes. The shockwave 
profile is 

$(r ,y = 0) ~ R s \j^^Cie~ Mrr ~ e~ Mir (2-10) 

where M\ is the mass of the lightest graviton excitation: if we are reducing gravity on the 
IR brane, gravity has a nonzero KK mass M\. Then one finds the gauge theory cross section 

ny 

<7gauge - ATr[—Zn[0.5v'sMiG/t]] 2 (2.11) 

where K is a numerical constant depending on the details of the gravity dual, M\ translates 
to the mass of the lightest QCD excitation and the constant multiplying ^/s in the log cannot 
be taken too seriously, as the subleading behaviour of a is modified anyway. 

The description of this last Froissart behaviour was shown in (3j to match exactly to the 
description of the saturation in Heisenberg’s model |SJ. Heisenberg says that at high enough 
energies the Lorentz-contracted hadrons colliding will look like shockwaves, characterized 
by a transverse size ~ 1 /M H . Moreover, at high enough energies, the hadrons effectively 
“dissolve” in the pion field, which also becomes Lorentz contracted to a shockwave, with 
transverse size characterized by l/m v . 

Heisenberg assumes that the “degree of inelasticity” a {=£/yfs= energy loss/collision 
energy) behaves like the overlap of pion wavefunctions, and the pion wavefunction behaves 
like ~ e~ m ” nr . Then one finds that the behaviour of the QCD cross section is given by 

o = vr b max (s) 2 ~ —jin? ^ (2.12) 

m- < E 0 > 

where < Eq > is the average emitted pion energy. But the average emitted pion energy is 
found to increase linearly with energy for a free pion action or for an action of the A0 4 type. 
Instead, Heisenberg finds that one needs a nonlinear action in order to get an approximately 
constant < Eq >. He takes the DBI action in (d and then he gets < Eq >~ m n . 

The same kind of picture appears in the gravity dual, where we collide particles, charac¬ 
terized by some size, but in the high energy limit we find that the only thing of relevance is 
the gravity field, and we effectively collide gravitational shockwaves, also characterized by a 
shockwave profile that has a characteristic size M 1; the scale of the lightest excitation. 

Until now, we talked in the gravity dual about pure gravity, corresponding to pure gauge 
theory. But in reality QCD has pions (Goldstone bosons of chiral symmetry) which are 
much lighter than the scale of the lightest glueball. The discussion of the gravity dual before 
the onset of the Froissart behaviour is the same, as it is governed still by gravity producing 
black holes. But if one has a Goldstone boson in the gauge theory, it can be modelled by 
a radion in the gravity dual (the position of the IR brane can be made dynamical). If the 


5 





radion mass is smaller than M\ (the KK graviton mass), the radion will dominate at large 
enough s (dominated by the IR). The scattering will then produce local brane bending, and 
the brane will bend significantly, entering the effective scattering region of the gravity dual 
before the black holes created in this region will become large enough, and thus the cross 
section at infinite s will be dominated by brane bending. If the radion is more massive than 
gravity {Mi), the brane bending will be small, and the created black holes will reach the 
brane before the brane reaches the scattering region, and the cross section at infinite s will 
still be dominated by black holes. 

In the physical case of lighter radion (thus lighter pion), at least heuristically the same 
Froissart behaviour will apply, as shown in m We do not have at our disposal the rigorous 
arguments of [2j anymore, since those were based on powerful general relativity theorems 
about horizons, and the collision of shockwaves in scalar field theory is still too hard to solve 
explicitly, but we will still rely on our gravity intuition for guidance. 

One observes that the action for the radion, considered as a brane moving in one (almost 
flat) direction, is in fact the DBI action 

S = T ( 3) J d A x^J 1 + (V0 2 /A 2 (2.13) 

where X is the radion (brane position) and A is related to the string scale. This is exactly 
the action taken by Heisenberg im at zero mass, but one could consider as a radion 
stabilization mechanism giving it a mass according to dm Moreover, if the IR brane is 
considered to be a D-brane probe, it will also have a U(l) gauge field on it. On static 
solutions (time independent) and at zero magnetic field, the action can be taken as 

S = T C3) J d A xiJl + (VW) 2 /A? - (V0) 2 /A| (2.14) 

where 0 is an electric potential and Ai and A 2 can be a priori different. We have to remember 
though that this is the D-brane action in a flat extra dimension, thus the curvature of the 
gravity dual space will induce higher order corrections to the DBI action. We will discuss 
them later. 

We are thus driven to study the DBI actions with both signs inside the square root, 
corresponding to either real scalar or electric potential (coming from the original Born-Infeld 
action). 

3 DBI actions and Blon solutions 

In this section we look at DBI actions and their Blon solutions and try to connect to the 
picture from the last section. 

We saw that both in the Heisenberg model and in the gravity dual description we are 
driven to consider DBI actions. In the gravity dual shockwave collisions produced black 
holes, and equivalently we proposed that in QCD pion field shockwave collisions should 
produce “fireballs”, solutions with horizons radiating thermally. The collision process is 
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quite complicated, as we saw in the gravity dual case, where we could prove a black hole 
forms, but we couldn’t calculate its metric, so we took as an approximation that we create a 
static spherically symmetric black hole. In reality however we have a dynamic complicated 
process. Similarly now we will take as an approximation that we create static spherically 
symmetric “fireballs,” thus in this section we will study such solutions. 

Moreover, one can ask also whether one can understand the nucleons as being modelled 
by solutions to the same DBI actions. We will thus look at all static spherically symmetric 
solutions of the DBI actions. In this section we will set the DBI scales A] and A 2 to 1 for 
simplicity. 

DBI action for 4d YM. The original action of Born and Infold m is 


h ' F » vFtlV - { F ^* FiLV )2 = J 1 _ E 2 + B 2 - (E ■ B) 2 


£= \l-det(r)^+ ~^=) = yi+ 2 


Here = Foil Fj = 6ijkBk. One defines in the usual way 

dC 


D = 


dE 


Then at B=0, the held equation is 


V -D = p 

as in electromagnetism (Maxwell), just that now we have a different definition for D 

E - D 

D = - , => E = 


(3.1) 

(3.2) 

(3.3) 

(3.4) 


Vi - ~e 2 VTTW 

Defining the electric potential 0 by E = — V0 we get the DBI-electric lagrangean 

£ = v^(Vf (3-5) 

4d scalar DBI. We follow the above analysis closely. In the static case (() t = 0), we 
define 

F = VX (3.6) 


The Lagrangean is then 


Defining 


£= v /-det(i,^+a„Xa„A')= \/l+(3„V) 2 = V! + f2 

c = sc 

the equation of motion is 


dF 
V -C = p 


(3.7) 

(3.8) 

(3.9) 
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as for the free theory! (analogous to the fact that for B1 we had the same equation as for 
the Maxwell case) except now the definition of C is different 


C = —=> F = —(3.10) 

Vl + F 2 V1 - C 2 

Blons and catenoids 

For a treatment of solutions to the scalar and electric DBI action see [221 ESI EH- The 
electric Blon (solution to electric DBI), originally found by Born and Infeld, is given by 


<j>{r) = C 

with the asymptotics (eqs 72-75 in 

C 


dx 


\JC 2 + x 4 


0(r) ~ —, r —> oo; 0(r) ~ const. — r, r 
r 


(3.11) 


(3.12) 


From the asymptotics at infinity, we see that C = g=electric charge. The Blon is found by 
noting that the solution to (D with delta function source is the same as for the Maxwell 
theory in terms of D, and then finding </>. As Born and Infeld noted, the advantage is now 
that the solution is nonsingular in terms of E, which reaches its maximum of E = 1 at r=0. 
Consequently, also the energy of this solution (with delta function source) is finite, unlike 
the case of Maxwell theory. 

Analogously, one finds the “catenoid” solution, i.e. solution to the scalar DBI action, as 


X(r) = C 



dx 

Vx 4 - C 2 


(3.13) 


But this has a horizon at r = r 0 = \[C, where F = V X diverges, even though X remains 
finite. 

Solutions to the D-brane action 

As shown in (221, one can embed the previous solutions into the U(l) D-brane action, 
that on static solutions (and at zero magnetic held) is 


S = T(3) f d'xsj (1 - (V<« 2 )(1 + (VX) ! ) + (V0 • VX)> 


(3.14) 


For a review of solutions to this action, see also ( 21123 - On top of those, the D-brane action 
also has BPS solutions, which were shown to correspond to (C — C = q) 

(j) = X = - (3.15) 

r 

and are understood in string theory as fundamental strings attached to the D brane and 
stretching all the way to infinity (BPS Blons). It has a singularity at r=0, where the 
fundamental string is attached, and has infinite energy, because the fundamental string has 
finite tension and infinite length. 











The Blon solution however has no singularity nor a horizon at r = \fC. One solves the 
equation for D in the usual way (D = X(C/r)), and at r = y/C, \D\ — 1 and \E\ — 1 /v/2, 
and there is no singularity: the electric potential <f> is continuous and its derivative is also. 
The Blon has a finite total energy, which is the reason Born and Infeld put it forward in 
the first place, as they wanted to have an electron (solution with delta function source) with 
finite classical self-energy. In string theory it is not clear what its interpretation is, but one 
might not worry about that too much, since it carries no topological charge, so one could 
maybe have doubts about its stability. The simplest possibility for an interpretation is of a 
string going through the D-brane, such as not to excite the scalar field. 

An interesting property of Blons is that like charges repel each other, and opposite 
charges attract (just as in Maxwell theory), and one is able to write down explicitly a Bionic 
crystal solution (first done by Hoppe [25], see also j2Sj)j in which each charge is surrounded 
by opposite charges, like in a NaCl crystal of Maxwell theory. 

For the “catenoid” solution, at r = \[C, X is finite but not its derivative, and we will 
call this surface a horizon. We will note later on that this name is justified, and it is very 
like the horizon of a black hole, but for the moment it just indicates the aparent singularity 
of the solution. Putting back momentarily an energy scale /3 2 , i.e. from 

S = /r 2 J [y/lTWW? - !] + f X(C5{r)) (3.16) 

we have 

C, = ft-; F t = , S (3.17) 

A - 

We can ask the question: can we continue the solution beyond this aparent singularity? 

Unlike for a black hole, it was argued in [22] that in this case the only analytic continuation 
that makes sense is to glue another solution with a different asymptotic space, creating 
an Einstein-Rosen bridge. In string theory, this corresponds to a brane- antibrane pair, 
connected by a fundamental string of length L, as in fig{T] 

But if we think of the catenoid as being a metastable solution, created in a collision of 
some type, as we will want later on, one should not have a full extra anti-D-brane, but one 
should at least be able to continue into something that looks locally as an extra brane (a 
brane “bubble”) as in figJ21 if not in a different way. 

For a black hole, when we reach the horizon, we can analytically continue in two ways: 
we can either continue to a singularity, OR glue another asymptotic region to create an 
Einstein-Rosen bridge. So if the catenoid is a scalar analog of a black hole, we would think 
we should be able to continue the solution to r=0 without making an Einstein-Rosen bridge. 

However if we naively continue F through r = \f(J using the same relation relating it to 
C (where C is the solution to V-C = C5(r)), it becomes complex, thus X becomes complex. 
So the only possibility of keeping F and X real is to change the sign of ( dX ) 2 in the action, 
effectively changing X into a (j) (electric potential)-type variable. Then the continuation is 
analytic in the sense that C has the same (analytic) expression inside and outside r = \/~C 
(for it, nothing interesting happens at r = VC). But it is not clear what would be the 
physical significance of changing the sign of ( dX ) 2 in the action at the horizon. 
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anti-D-brane 




Figure 1: D-brane anti-D-brane system connected by a string: two catenoid solutions with 
different asymptotic regions connected on the horizon 



Figure 2: A single D-brane, the continuation past the horizon can be a brane “bubble” 
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Note that 


0(VC) = C 


dx 


<Vc Vx 4 + C 2 
thus 0(0) = 20(v / C'), whereas 

X(VC) ^ r dx 


= V c2_Ewr = c 


~Vc 


dx 


71 


\J x 4 + C 2 


= C 


'Vc V‘- 


x* 


w = 


c=c 


(3.18) 


(3.19) 


so gluing the X solution to a 0-type solution inside the horizon would give a finite displace¬ 
ment X(0) at r=0 (A" is the position of the D-brane in the transverse direction). 


For the Blon action 


(V0) 2 , we have also the solution 


0(r) = C 


dx 


VC 2 


x^ 


(3.20) 


which is just the continuation of the “catenoid” solution to r < ro, and redefining X —> iX, 
thus effectively changing the action. This solution also has a horizon at r = ro, where 0 is 
finite, but 0' is infinite and negative. However, this solution has imaginary action (or energy) 
(the square root is negative), but is a real slution to the equations of motion. 

Finally, note that the Blon and catenoid solutions are defined up to a sign (±), and an 
additive constant. One could construct more solution that are not real. For instance, ix the 
Blon solution is a solution to the scalar (“catenoid”) action, but that is a trivial observation, 
reflecting the fact that redefining the scalar by an i takes us from one action to the other. 

For a general static solution of the D-brane action studied in j2|, the equations of motion 
can be written as 


V ■ g p n = V ■ 


E( 1 + (VX) 2 ) - VX(E ■ VX) 


= 0 


'(1 - £ 2 )( 1 + (VX) 2 ) + (E • VX) 2 

vx + 3 2 n(n • vx) 


1 + (VX) 2 + g 2 U 2 + g 2 (Jl ■ VX) 2 
VX(1 - E 2 ) + E(E ■ VX) 

(1 - E 2 )( 1 + (VX) 2 ) + (E ■ VX) 2 


= 0 or 


= 0 


(3.21) 


On a spherically symmetric solution (with E = — V0), the equations become 
r r 2 0' „ „ r r 2 X' 


' y/1 - 0 /2 + X' 2 


r = o; 


' y/l ~ 0' 2 + X' 2 


r = o 


(3.22) 


The solutions of these equations are obtained by integrating them with constants C and 
C, respectively: 


./2 i + x' 2 x , 2= l -4>' 2 


1 + r 4 /C 2 ’ 


r 4 /C 2 - 1 


(3.23) 
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from which we obtain the general solution 


A' = 


C 

Vr 4 + C 2 -C 2 
C 

Vr 4 + C 2 - C 2 


(3.24) 


Thus we see that for C > C, both 0 and X look like the Blon, being everywhere defined and 
having finite derivatives at zero. For C < C, both 0 and X look like the catenoid, having a 
horizon at a finite r, equal to ( C 2 — C* 2 ) 1//4 . At C = C = q we have the BPS Blon, which 
blows up at l—O. 

At r=0, we can specify 0(0) = 0'(r = 0) and A(0) = X'[r = 0), but we have to satisfy 


(1 - 0(O) 2 + A(O) 2 )0(O) = (1 - 0(O) 2 + A(0) 2 )A(0) = 0 (3.25) 

i.e. either 0(0) = A(0) = 0, or 0(O) 2 — A(0) 2 = 1 (for the BPS Blon 0(0) = A(0) = oo, so 
that solves it). We can check explicitly from the general solution that this equation is true. 
The Blon has 0(0) = 1, thus we can think of any solution with X as a perturbation around 
the Blon. Note that this means that the square root in the action is equal to zero at r = 0. 
Only for the BPS Blon the square root at r = 0 is equal to 1. 

At r = oo, by putting 

0 ~ — X ~ — (3.26) 

we find from the equations of motion that n = m = 2, as for the separate 0 and X theories, 
thus the interaction doesn’t change that, as it should. 

As we saw, if we put 0 = 0, any solution will have a horizon (the unique solution is of 
catenoid type, which has a horizon). We check that indeed, we can’t have a real solution 
defined at r=0, since then — A(0) 2 = 1. However, if 0 ^ 0, we see that the solution has a 
finite A(0) = y/0(O) 2 — 1, and the solution can be extended to infinity without encountering 
a horizon. This is due to the fact that we can think of these solutions as adding more electric 
field at r=0 (0(0) — 1) to the Blon, in order to compensate for the added scalar (A"(0)). 
The Hamiltionian of a static configuration (A = 0) is 

1 + (VA) 2 

(1 - E 2 )( 1 + (VA) 2 ) + (E • VA) 2 


1] (3.27) 



For the general spherically symmetric 4d static configurations above it is 


„ 1 f ,, r 1 + A' 2 

E — - I d 3 x [— 

9 J y/l - 0 /2 + A' 2 


n 1 /‘ ,3 r rA + ° 2 

- 1 = - / o ,- == 

9.J Vr 4 (r 4 + C 2 - C 2 ) 


- 1 ] 


(3.28) 


Specifically, for the catenoid (purely scalar), 


gE x = I r 2 dr(V 1 + A' 2 - 1) = C 3/2 I 


(3.29) 


'r o 
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where 


/ 



1] ~ 0.770343 


(3.30) 


Note that then r 0 = \/~C oc (. Ex ) 1 ^ 3 , same as for an object of constant energy density. 
However, the energy density diverges near the horizon, showing that a significant fraction of 
the energy is concentrated near it. For the Blon we have 


9E4, 



r 2 dr[ 


1 

\/l - 0 /2 



dx 

x 2 + a/x 4 + 1 


(r[i/4]) 2 C 3 /2 

6v^r 


(3.31) 


and now the energy density diverges at r=0, thus a significant portion of the energy is 
situated near the origin. Then for the general solution with horizon {C > C ) we have 


gE 


(~i2 

r 3 0 I+~V« 
r 0 


£ 15 / 4 ] , 

r[3/4]’ 


r o = 


G 2 


(3.32) 


and for the general solution with no horizon [C > C ) 


gE 


(r[1/4 ]) 2 fg G 2 

2y/irf Q L 3 2 


f 0 4 ee G 2 


C 2 


(3.33) 


For the Blon C is an asymptotic U(l) charge, thus will be quantized in the quantum theory, 
and for the catenoid G is an asymptotic scalar charge, thus again we expect it to be quantized. 
We see that both the catenoid and the Blon energies go like charge to the power 3/2, thus 
higher charge objects will be unstable towards decay onto the lower charge ones (E(Q 1 + 
Q 2 ) > E(Q 1 ) + E(Q 2 )) This is also valid for the general solution with both C and C nonzero 
(if C oc C oc n, then E oc n 3//2 always). 

Due to the nonlinearity of the action, it is hard to construct explicit (separated) multi 
center solutions. The classical interaction potential between two solutions would be then 
given by the difference between the energy of the two-center solutions and the individual 
energies of the single center solutions, 


E(R) = E( 1, r 0 = 0; 2, r 0 = R)-E 1 ~ E 2 


(3.34) 


Given that we can’t construct the multicenter solutions, we can’t calculate the form of the 
potential, but we can say something about the asymptotic features. At large distances, 
the interaction becomes (free) Maxwell electromagnetism plus free scalar, thus the potential 
between two Blons will be 

E^R) ~ (3.35) 

XL 

where Qi and Q 2 are the electric charges {C\ and C 2 ), thus repulsive if Q\Q 2 > 0 and 
attractive if Q 1 Q 2 < 0. And the potential between two catenoids will be 

E X {R) - (3.36) 
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thus attractive for like scalar charges Q\ and Q 2 (£ 1 , 62 ). For a general solution, it should 
be the sum of the two potentials. 

However, at small enough distances, the two center solution will look approximately like 
a single center solution with charges equal to the sum of the individual charges, thus the 
potential will always be repulsive! That is true both for the Blons, and for the catenoids, 
thus for the general solution, since in all cases the energy goes like charge to the power 3/2. 
This is of course for the case when the interacting objects are of the same type (same ratio 
C/C). 

The only exception is the BPS Blon, for which the energy of two center BPS objects 
exactly equals the sum of the individual BPS objects, i.e. the potential is zero. In that case, 
the potential is zero at infinity (E^ cancels against E x ), but what about at zero? 

The energy of the BPS Blon is 

gE= I d 3 x(VX) 2 = jf r 2 dr[ r * + ^ 2 _ ^ = C 2 ^ (3.37) 

thus divergent, and the divergence is the same as in the free scalar theory (or free electro¬ 
magnetism). Because the energy is divergent, we can’t draw any conclusions from the fact 
that it behaves like the square of the charge. I 11 fact, the energy of the multicenter solution 
is exactly equal to the sum of the individual solutions, due to the BPS property. The point 
is here that one has to regularize the infinite energy, and if one introduces a lower cut-off for 
r, r > 5, while keeping also the value of the scalar field, X(5), fixed, one gets 

C 2 

gE ~ — oc CX(S) (3.38) 

0 

thus energy that is linear in the charge C. The divergence signals the fact that we have 
to take into account quantum theory, and at the quantum level we know the BPS Blon is 
a BPS string, thus stable, and in string theory X(S) has physical significance. Thus the 
asumption is that 77(5) is kept fixed, i.e. Xt(<5) = X 2 (5) = X 1 + 2 (5), and one just adds the 
coefficient C. For a catenoid we have also X(f 0 ) = <77 1 // 2 -y/^rF[5/4]/T[3/4] ~ 1.31103C' 1 / 2 , thus 
E oc CX(r 0 ), but now the energy is finite, and we can treat this object classically, and there 
is no reason to keep X (fo) fixed, but we can let it vary with f 0 , and there is no interpretation 
in terms of a string (fo is just a finite “thickness” of the solution). 

The same situation of diverging energy would happen if we tried to apply the above 
logic to a free theory like Maxwell electromagnetism or free scalar theory. For example, for 
electromagnetism, the classical potential would be 


E(R) = / d 3 x[E 2 ( 1, r = 0; 2, r = R) - El - E\ 


= / r 


2 . U QC . Qi{f - K)/\r - R\ 2 Qrr 2 


J dr[(- 


+ 


Rf 


L )-(^r-( 


Q 2 (R-R)/\r-R \, 2 


Rf 


) 2 ] 


= 2 Q 1 Q 2 r dr 


2 f ■ (f— R)/\f— R\ 


r 2 (r - Rf 


(3.39) 
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and this diverges, and we can only draw the conclusion that E(R) ~ Q1Q2/R by scaling, 
which is correct, but nothing about the coefficient or its sign. In fact exactly the same 
calculation applies for the free scalar, but for scalars the true sign is opposite! 

4 Boosting Blons; comparison with gravity 

In order to have a consistent effective field theory picture for high energy scattering, we 
should boost the solutions corresponding to the nucleons, creating shockwaves in the pion 
held. The collision of such shockwaves should create “fireball”-like objects. In this section we 
will thus analyze the boosting of the DBI solutions from last section. We will then compare 
with the gravity dual case. 

Boosted “catenoid” 

If one boosts a scalar static solution f(x,y,z ) , by definition the boosted solution is 
g{x',y ', z!,t') = f(x,y,z), thus for a spherically symmetric solution f(\jx 2 + y 2 + z 2 ) we 


have 



and it thus satisfies the equations 

dz'9 = ld z f\ z=l{z ,_ vt) ] d y rg = d y f\ ; d x ,g = d x f\\ d t >g = -'yvdj] (4.2) 

Let’s then boost the “catenoid”, 



(4.3) 


This is a true scalar, so we just replace r with 



(4.4) 


Then after the boost 

ax c f ox c 


7 2 (z — vt) 


dr ~ yV 4 — /3 2 C 2 r 1 dz ~ ^ r 4 - f3 2 C 2 
dX _ _ C 7 2 v{z~vt) 
dt ^r 4 - (5 2 C 2 r 



r 4 - p 2 C 2 


r 


as before boost (4.5) 


Thus in the limit r->lwe have the shockwave-like solution 


X = X(r), z — t 
0; otherwise 

(one could write X = X(r)5(z — f)/<5(0)). But at finite v, we have 
dX/dz 7 (z- vt) 

, r = 7 -t-> 7 —> oc (it as v —> 1 , 7 (z — vt) ^ r 


ij” 'f 


(4.6) 


(4.7) 
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Let’s check the Lorentz invariance of the held distribution by looking at its energy, and 
proving E = ^E 0 . From 

S = — J d 3 x dt(^l + (d^X) 2 - 1 ) (4.8) 

we hnd 

1 + {XX) 2 

^i + {xx) 2 -x 2 

Then before the boost, 


i] 


(4.9) 



Eq — / 47ir 2 dr[\ 1 + 


C 2 


/ 3 2 C 2 


— 1] = 2nrdrdz[\ 1 + 


C 2 


/3 2 C 2 


- 1 ] 


(4.10) 


and after the boost (changing the variable of integration from z to z' = 7 (z — vt )), we get 

c 2 


E — — I 27T fdrdz' [ t /1 + 


C 2 


1 + 


r 4 —/ 3 2 C 2 


r 2 +7 2 z' 2 


7 


r 4 - p 2 C 2 1 + 


c 2 

r 4 -/3 2 <72 


1 ] 


(4.11) 


where r 2 = f 2 + £ /2 . We can see at most that at large 7 , E goes like 7 2 /7 = 7 , as it should. 

Boosted Blon 

Now the scalar is actually the electric held, the zeroth component of a vector, thus the 
action is 


S = — d A x dt[\/1 - (V(f> - A) 2 + (V x A) 2 - (C Vcj) - A) • (V x A)) 2 - 1] (4.12) 


and becomes 


S 


d 3 x dt[\J 1 — ( dcj )) 2 — 1 ] 


(4.13) 


only when A = 0, 0 = 0, thus we have to boost differently, taking into account that (j) is not 
a scalar anymore, but now (0, A) is a vector. 

Then on a static, purely electric configuration like the Blon {F 0i = E i: Fij = e^Bk) we 
get 


E' z (x', y, z\ t') = E z (x, y, z); E' x (x' M ) = jE x (x^); E' y {x' y ) = 7 E y (x^ 

ByW = 7 vE x (Xf,); B' x (x'J = -7 vE y (x^)] B' z = 0 (4.14) 

Taking Ei = —<9i/(\Ac 2 + y 2 + z 2 ) and boosting it to obtain E[, = d f >A), — di'g(x',y, z',t'), 
and if we also boost (/, A) as a vector (i.e. giving an extra condition, that we don’t change 
gauge when boosting, which is a nontrivial condition, as this is not the case for gravity, see 
below), we obtain 

g(x, y , z, t) = A' 0 = 7 (A) + vA z ) = jA 0 = 7 f{\jx 2 + y 2 + ^ ■) (4.15) 
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Specifically, for the Blon we have 


%= , = 7(2 - Pt) ; % = -=£==£■, B' t = - C (4.16) 

Vr 4 + /3 2 C 2 r- -/r 4 + /^C 2 r y ^r 4 + /3 2 C 2 r 

thus in the large 7 limit, 

Sr = / r , • 7 -> 00 if liz - vt) ~ f; B' g ^ E' r (4.17) 

71 - - uf) 

If we boost </> as the 0 component of a vector, at n —> 1 we get that ft = 0 if z 7 ^ t, but now 
4>(z = t) ~ 7 —> 00 . Moreover, the width in z — vt is ~ I/7, thus 

4>'(f, z, t ) = a(j)(r)8(z — t ) (4-18) 


with a a number. 

The extremal Blon of Callan and Maldacena 223 is a solution to the action in m . 
The static BPS solution is 

X = i ■, E = VX (4.19) 

| r — ro\ p z 

thus both X and E are singular at r=0. Boosting this will generate the same X as for the 
catenoid (with horizon at r= 0 ), but now we will have both electric and magnetic fields as 
well (thus the action in (18.141) is not valid anymore, it needs to be “boosted” as well, i.e. the 
B dependence specified). 

Note that when we boost the black hole solution to obtain the A-S solution, we write 

2MG 

9tt(xy', z', t')) = 1- ; g rr = 1 /g tt (4.20) 

A' 2 + y' 2 + 

and then we transform gtt,9rr to g' )v , (and take the limit v —> 1 ), which would be the 
equivalent of calculating only E', not g. Moreover, one then has to transform the coordinates 
also to reach the system where the A-S looks simple. It could be that the same is required 
here (for Blons and catenoids), but we will not pursue this further. 

5 Blon solutions and “dumb holes” as fireballs 

In |ll)j it was found that there are analogs of black holes in hydrodynamics, dubbed “dumb 
holes”. When a configuration of fluid moves at ultrasonic speed it creates horizons that 
radiate thermally, analogous to the thermal horizons of black holes. The hydrodynamics 
equations are written for a potential flow, in terms of a scalar potential $. On the other 
hand, we want to obtain also a thermal horizons in the collision of scalar (<h) shockwaves in 
DBI theory. We will see that in fact the equations in the two cases are completely similar. 
We will thus first review the “dumb holes” and then completely parallel the calculation 
for the DBI scalar, calculating the temperature of scalar “fireballs”. For the “fireballs” we 
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will take the static spherically symmetric solutions studied in section 3, and we will discuss 
towards the end whether this is a good approximation to the real dynamical situation. 

For ultrasonic fluid flow, the surface where v = c (the velocity of particles reaches the 
sound velocity) is a horizon that radiates particles thermally, and can be mapped to a black 
hole. The fluid has an equation of state p = p(p) and the fluid motion is irrotational 
V x v = 0, such that one can write v = V<h (potential flow). The speed of sound is defined 
by c 2 = dp/dp. 

The equations of motion are the local pressure equation and the continuity equation, i.e. 


p[y + [y • V)u) = -Vp(p); p + V ■ (pv) = 0 
The first equation is integrated to the Bernoulli equation 

* + ^ + ft(p) = 0; 

2 J p 


(5.1) 


(5.2) 


If we now calculate the fluctuation equations, in variables 0 = <5$ and 0 = 5p/p, we get 


0 + v ■ V0 + C 2 0 = 0 

i^r + v- V)0 + (V0) ■ V In p + V 2 0 = 0 
dt 


Eliminating 0 we get 


1(0 + v- V + (V ■ *0)4(4 + »' %> - Iv(pV^) = 0 

p dt c z dt p 


(5.3) 


(5.4) 


which is exactly the equation of motion for a scalar field in a curved spacetime, d^y/gg^d^cj) = 
0 if the metric is given by 


y/agT = p 


J_ v*_ 

c 2 . c 2 

yi v' L vi £ij 


(5.5) 


which implies in 4d (after finding g lw and defining a new time coordinate by dr = dt + 
v l dx */(c 2 — v 2 )) 


2 > 2 

or, in the case of a radial flow 


ds 1 = ^((c 2 - v 2 )dr 2 - (S‘ j + 0^4 
c c z — v z 


)dx l dx j ] 


i 2 P 17-1 2 / 2\ 2 j 2 ^V 

ds = -[(1 — v /c )c dr — 


C" 1 — v 2 / c 2 

In the new coordinates, the scalar wave equation is 


- r W] 


^0 7 —^ 2 / 2 ^°^ 

1 - v z /c z 


v l v J 


S ij )d j( j) = 0 


(5.6) 


(5.7) 


(5.8) 


At constant (or neglijible variation of) p, the fluctuation equation in the original coordi¬ 
nates can be written as 


[d t + (V$) ■ V + V 2 <I>Pf + (V<h) ■ V]<5$ - c 2 'V 2 8$ = 0 


(5.9) 
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The horizon is the surface where v = c ((V<I>) 2 = c 2 ), thus the metric becomes singular 
and then in the above equation we have 

[d 2 t + ((V$) 2 - c 2 )V 2 + = 0 (5.10) 

thus the coefficient of k 2 in a fluctuation mode k changes sign. 

For static fluctuations (d t = 0), the equation is 

—c 2 V 2 h<f> + (V 2 $)(Vd>) ■ V(hh + (V$) • V((Vd>) • W<f>) = 0 (5.11) 


DBI Lagrangean 


The DBI scalar has a Lagrangean hl + (V<f>) 2 — (d t $) 2 . Here we denote the scalar by 
just to underline the analogy with the hydrodynamics case. It gives the fluctuation equation 
on static solutions (<9yl> = 0, d t 5$ ^ 0): 


<5[V 


V<f> 


1 + (V4>) 2 
1 


1 + (V4>) 2 


=< 9 2 < 5 $ 


-V- 


1 + (Vd>) 2 
(1 + (Vd>) 2 ) 3 / 2 


df + V ■ 


Vh<I> 


0 - 


1 + (Vd>) 2 
((Vd>) • W<F)] = 0 


(5.12) 


We see that now the equation is very similar to the one above. The horizon is again 
where the coefficient of V 2 <5<I> changes sign, which now means 


(V4f 


d+(v$m= 


(5.13) 


whose solution is |V<L| = oo, which is what we called a horizon for the “catenoid” anyway! 
Before we thought of the |V$| = oo solution as a horizon just because of the singularity, 
but now we see that it is the exact analog of the dumb hole horizon, thus of the black hole 
horizon. 

In fact, we can do more. The fluctuation equation m can be rewritten in a suggestive 
form as (using that <%$ = 0) 


-d t 


1 + ( V $) 2 


--d t 5$ + di 


1 + (V<F) 2 


--{S 






d J i> 


z)dj5$ = 0 (5.14) 


1 + (V4>) 2 Jl + (Vd>) 


which is the same as the black hole equation in the new coordinates (EH) with the ident.id¬ 
eation 


c 2 = 1 + (Vd>) 2 , p = 


1 + (V$) 2 


v 

c 


d l (l> 


1 + (Vd>) 2 


(he., v l = d *$) 


(5.15) 
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Of course, just because v=c when | V<F| = oo it. is not a good reason for it to be a thermal 
horizon. There are many hydrodynamic surfaces where v = c that do not radiate thermally 
as “dumb holes”, just as there are horizons of finite area, for extremal black holes, that 
nevertheless are at T=0. The key point is that the variation at the horizon, dv/dr\ v=c for 
dumb holes and d r gtt\ r =r H oc k for black holes that gives the temperature of the horizon. But 
the point is that neither for the black holes or the dumb holes, the presence of a singularity is 
not required, the thermality is a property of the horizon itself. What lies behind the horizon 
is of no consequence. 

Information trapping and horizons 

One should observe though that the horizon of the scalar fireball doesn’t have the black 
hole horizon property that nothing can escape from the horizon. It has that property only 
for scalar held excitations and only in a very limited sense. In general, excitations cannot 
come out of the horizon because of the same reason that that the horizon radiates thermally 
(the coefficient of V 2 </>, proportional to v 2 — c 2 , becomes zero at the horizon in a specific 
way). This is seen as follows. 

For a Schwarzschild black hole, the infinite time delay necessary for a massless particle 
to come out, from the perspective of an outside observer is due to equating ds 2 = 0 and 
finding that then 



rR 


dr 


c(l — 2 M/r) 


oo as R —> r H 


(5.16) 


For the metric in m this becomes f dr/[c(l — v 2 /c 2 )\ being divergent at the horizon, which 
is easily seen to be due to the fact that d[c(l — v 2 /c 2 )\/dr is a nonzero constant at the horizon. 
This is the same condition that we will fold below for the nonzero horizon temperature, if p 
is a constant or satisfies pc( 1 — v 2 /c 2 ) ^ constant (in our case p is not constant but satisfies 
the latter property). Thus the infinite time delay and finite nonzero temperature are almost 
always related, as expected, and if it’s true, scalar excitations cannot get out. However, in 
our case, d[c( 1 — v 2 /c 2 )\/dr is infinite at the horizon, implying a finite time delay, and we 
will fold below an infinite temperature as well. 

But any other (nonscalar) excitation can always come out. We might ask how come <F 
fluctuations can’t get out if d[c( 1 — v 2 /c 2 )]/dr is nonzero and finite at the horizon, after all, 
we are in flat space so it takes a finite time for radiation to come from the horizon. But 
the point is that the scalar fluctuation equation is the same as for a black hole, meaning 
the characteristic equation for a scalar held perturbation will propagate as a function of the 
asymptotic coordinates r, t the same way as from a black hole. In other words, the phase 
and group velocities ( c p h and c gr ) of scalar excitations would tend to zero in those cases and 
information could not be exchanged with the inside of the horizon. 

In our case, substituting = A exp(i(cot — kr )) (spherical waves) in the perturbation 
equation dSH2D one obtains 


oj 2 = 


1 + <F /2 


{k 2 - 3 k 


l + $' 2 


) = k 


4 4 

2 r ~ r 0 


+ 6 k^r 


(5.17) 


thus we see that the coefficient of k 2 becomes zero at the horizon and the coefficient of k is 
finite. One can easily check that if for our solution we had ~ l/(r — r 0 ) near the horizon 
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(which would imply both finite nonzero temperature and infinite “geodesic” time delay) we 
would get also the coefficient of k to go to zero at the horizon, but slower than the coefficient 
of k 2 . As it is, we get the phase and group velocities 





6 

-— as r -> r 0 
kr o 


duj a + b/(2k) 1 l~b 1 / 6 

C « r = lk = a + b/k ^ 2\l k ^ 2\J hT aSI ^ r » 


(5.18) 


and we see that infinite k modes have zero phase and group velocities. This is not a relativistic 
formula (u = ck, c= constant), but not a nonrelativistic one either (which would generically 
be of the type c = cq + ak + bk 2 + ...), but comes from the gravitational (black hole) 
background. The fact that one could apparently have c p h > 1 and c gr > 1 is an artefact of 
our approximation of spherical waves, which clearly only works at the horizon if k » l/r 0 
(A <C r 0 ). 

Thus high energy perturbations have phase and group velocities going to zero as 1 /Vk, 
hence for them the horizon is indeed impenetrable, and information can be exchanged with 
increasing difficulty with the inside of the horizon. If we would have infinite “geodesic” time 
delay ($' ~ l/(r — r 0 ) near the horizon), then the phase and group velocities would be zero 
at the horizon, thus all scalar information would be stuck at the horizon, exactly as for a 
black hole. Then the phase and group velocities would be again proportional to 1 /Vk, but 
would be multiplied by a factor that vanishes at the horizon. 

Finally, notice that if we take the original electromagnetic DBI action for 0 (for the Blon, 
with the minus sign in the square root), the sign of the two terms in (15.121) is the same, thus 
one can never cancel them against each other. Thus the original DBI action never admits 
horizons! 

Temperature calculation 

For a black hole in flat space, the horizon temperature is 

r = A (5-M) 

where k is the “surface gravity” of the horizon. For a static, spherically symmetric solution 
with only g rr (r ) and gtt(r) nontrivial (and possibly an r-dependent conformal factor for the 
sphere metric), one can easily calculate that 

(2k) 2 = lirn —(d r g tt ) 2 (5.20) 

horizon 

For a Schwarzschild black hole, g rr = g tt and we calculate that k = 1/(4 MG), as known. 
For the “dumb hole”, using the above map to a curved spacetime m we get that m 

(2k) 2 = {-d r (pc( 1 - v 2 /c 2 )]}\ =c =>T = ^-d r [pc( 1 - v 2 /c 2 )]\ v=c (5.21) 

p 4tt p 

Note that in Unruh’s case, where p and c are nonzero and finite at the horizon, one gets 
T = (dv/dr\ v=c )/(2 tt), but in our case that is not true. 
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For the catenoid, we have 


<£>' = 


D 


C/r 2 


\/l — D 2 71 - (C/r 2 ) 2 
and as we saw we map to the “dumb hole” calculation by defining the velocities (EED 

(V<F ) 2 


(5.22) 


v 2 /c 2 = [ 


(1 + (V $) 2 


-]; c 2 = 1 + <F /2 ; p = 1/Vl + <h /2 


(5.23) 


and then 


d r 1 


=* T = 


dr 1 + <h /2 
r 0 1 
4 (r - r 0 ) vrr 0 


2*: = |yrT^^[ T n_]| r „ o = |yrppA[F)u 0 = 


rp 4 
4(r - r 0 ) r 0 


(5.24) 


thus the temperature is infinite! 

Obviously, the temperature of a system cannot be infinite. There are several possible 
interpretations. First, we notice that the infinite square root factor in front of the temper¬ 
ature comes from y/1 + <f> /2 , which is equal to c using the map to black holes, and if c was 
a constant, we could have rescaled the time coordinate to get rid of it in the wave equa¬ 
tion. As it is, we cannot rescale it away. But \/l + <f> /2 is also the energy density of the 
solution, which cannot become infinite, so there will either be higher order corrections to 
the action preventing that, or else for the real dynamical (time dependent) process of high 
energy collisions one cannot approximate by the production of a static spherically symmetric 
solution. Thus at the classical level, the temperature could be made finite either by higher 
order corrections in the action, or by deviations from time independence and sphericity. As 
an example, the fluctuation equation around a time dependent solution is 


dr 


+V* 


+d t 


1 ■ <9*$ 

(a tJ - 


cP ( l> 




1 + (V<f >) 2 - <f > 2 yj 1 + (V $) 2 - <f > 2 \/ 1 + (V< 1>) 2 - $ 2 
1 

-d f 5$ 


1 _|_ ( v $) 2 - < f > 21 + ( v $) 2 - $ 2 




1 + (V<f >) 2 - $ 21 + (V 4>) 2 - $ 2 




-d t 


=(1 


$2 


1 _|_ (v $) 2 - $ 2 1 + (V $) 2 - <h 2 


^)d t s$ = 0 


(5.25) 


which will have a modified map to the black hole fluctuation equation. Finally, besides all this 
possible classical resolutions, an infinite temperature implies an infinite particle production, 
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thus quantum mechanics could also regulate this behaviour. In any case, we will assume 
that the infinite factor y/1 + <f>' 2 which is proportional, as noted, to the energy density of the 
solution, is regulated to some large but finite value, and calculate the mass scaling of the 
temperature coming from Tq. 

Restoring the energy scale Mp, we get 


1 1 Mp 3 

— = - oc -—— 

r 0 VC MV 3 


(5.26) 


(VC is related to the mass of the catenoid by M = IC 3 ^ 2 /g, as we saw in section 3.) 

This calculation is valid as long as ro < m^, thus for M < Mp/m After that we have 
to take rritj, into account. Similarly we then get 


T 



+ $' 2 + m 2 $ 2 


d 1 + m 2< h 2 
dr 1 + <f> /2 + m 2< h 2 


r=r i 


(5.27) 


Near the horizon, the equation of motion for the massive DBI scalar becomes 


r$"(l + m 2< f> 2 ) + 2$ /3 = 0 


(5.28) 


giving the approximate solution 


$(r ~ ri) ~ $1 + 


(5.29) 


This gives then 


T ~ 



r\ 


4(r — 7"i) 7T T\ 


(5.30) 


At large distances, $ ~ Ce~ mr /r, with C being some power law of the mass. Thus 
perturbatively, the horizon will be where <3> becomes of order 1, and if we would apply this 
perturbative formula we would get r\ ~ 1/m ln(C'm). Considering again that the infinite 
factor a/1 + < f )/2 + ^r 2 = \/l + m 2 $ 2 ^/r 1 /(4(r — rq)), proportional to the energy density at 
the horizon, will be regulated to a finite value and calculating only the mass scaling with rq 
we obtain 


T oc 


m 


ln(C'm) 


oc m 


(5.31) 


as we argued in |[6i we should get. 

Finally, if the temperature does become finite as we argued, by the arguments we already 
gave, light will take an infinite time to reach the horizon, a fact observed from geodesic 
arguments in the black hole background and from the vanishing of c p h and c gr in the scalar 
background. 


6 SU(2) DBI actions and Pavlovskii’s topological soli- 

ton 

Based purely on phenomenological grounds, Pavlovskii looked at a DBI action for the 
SU(2) pions of QCD. It contains an infinite number of higher derivative corrections to the 
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nonlinear sigma model action, that have a chance of having a Skyrme-like topological soli- 
ton. Such a soliton was indeed found. Here we will extend that discussion and look for 
modifications of the SU(2) actions that can fit our purposes. We will carefully look at the 
numerical solutions in order to generalize in the next section to the case of interest. 

The action E! writes is 


C = - 1 ] 


( 6 . 1 ) 


where = U 1 d IJ U, U = exp(i0 7r r// 7r ), and chooses a spherically symmetric configuration 
(“hedgehog”) 


V = e iF,r) « !f > ; ft = - 


( 6 . 2 ) 


One should note that the sign inside of the square root is different in El due to a local 
sign mistake. In the spherically symmetric solution used afterwards he has the same signs in 
terms of F, which is correct but because of the i’s in U, the sign changes between the above 
action and the action in terms of F (i.e. LjLj = — F' 2 + ...). 

Then 

U = = cos F(r) + in ■ crsin F(r) (6.3) 

and then for such a static configuration 


r sin F cos F sin F cos F sin 2 F 

Li = U^OiU = in ill ■ -) + %Oi - i - nje^kCTk 


L' z = - F /2 - 2 


sin 2 F 


(6.4) 


and the energy functional is (E — — S for static solutions) 


1 . F' 2 sin 2 F, 


E = Snf^l (1 — R)r z dr\ R= ,/!--(—+ ^ 


And the equation of motion is 

Explicitly, it gives 


j F 1 , sin 2 F 


R 


R 


(6.5) 


( 6 . 6 ) 


r 2 —— sin 2 F)F" + (2 rF' — sin 2 F) 


(3 2 


3 „/ 


—— (rF' 6 - F' 2 sin 2 F + -F' sirr F -- sin 2 F sin 2 F) — 0 

p 2 r r 2 


(6.7) 


For a good solution, F should go to zero at infinite r, thus the large r expansion of such 
a solution is found to be 

_, . d b c . __ 

F(r} = - 1 -—- 1 -— ( 6 - 8 ) 

V ' rpTl y>n-\-m y*n-\-p v y 
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and then from the equations of motion 


x Uj 

F ( r ) - I?- 


r 2 21 r 6 3 /3 2 r 8 

The first 2 terms in F(r) appear just from the equation r 2 F" + 2 rF' — sin2F 
would appear also from the “linear” action (from the SU(2) sigma model) 


5 = 


\dFf 


+ 


Sill 


2 F 


= - Tr- 


I 2 


(6.9) 
0, i.e. 

( 6 . 10 ) 


The first term appears when we approximate the sin to first order, and the second when we 
approximate to second order. 

By comparison, for the single scalar DBI action, the first two terms in the expansion 
would come from the equation r 2 F" + 2 rF' — rF' 3 / (5 2 , with the solution 


a a 3 
r 20 /3 2 r 5 


( 6 . 11 ) 


Notice that in first order, 


L„L„ — -id ■ d r dr) 2 = d^id^i 


( 6 . 12 ) 


which gives on a spherically symmetric solution (the above ansatz) 

F 2 

—F' 2 - 2— (6.13) 

Thus for three perturbative scalars, with action 

S = ^ J d 4 xd IJ n'd IJ n l (6.14) 

if one chooses a spherically symmetric (“hedgehog”) solution, with 7Tj = riiF(r),n = r/r, 
one gets 

S=\f d“x(F' 2 + 2^) (6.15) 

thus the difference in behaviour at infinity between the SU(2) case and the single scalar case 
comes not from the nonlinearity of the action (it couldn’t), but rather from the spherically 
symmetric ansatz. 

The above action has solutions that have T’(O) = Nn. Such a solution with N > 1, if it 
also goes to zero at infinity, will have a nonzero topological charge 

B = 2^ / d * xe i* Tr (LiL j L k ) (6.16) 

If T’(O) = Nn, then we can find from the equations of motion that F\ 0) = a is not 
constrained, thus parametrizing a set of solutions with F = Nn + ar + o(r 3 ). However, for 
a single value of a (negative) is the solution going to 0 at infinity, thus being topological. 
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In 103, the topological solution was found numerically, by imposing the good behaviour at 
infinity. For most solutions, one encounters a point where the solution becomes numerically 
unstable, at a nonzero r = r$ satisfying F(r$) = ±arcsin(/3ro) , where we can check that 
F'(ro) = 0. Physically, it is clear that thes solutions will be continued to F(0) = 0, F'( 0) > 0, 
thus carrying no topological charge. 

We have also analyzed numerically what happens when the asymptotics at zero are 
correct, i.e. F( 0) = tt (the other N’s are similar). Numerically, by trial and error, (everything 
is in units of (5) we find that the topological solution has a ~ —0.809466. Specifically, 
for F(0.01) = 7T — 0.00809466, F'(0.01) = —0.809466. Then F(300) ~ 0. Specifically, 
F(10) = 0.0706502 = b/r 2 => b = 7.065 and then F'(10) = -0.014102 ~ -26/r 3 , but the 
relation is not true anymore at r=100, showing that this is an approximate solution. This 
topological solution is given in figH3 with a detail in fig® 

Note that this value for a is very close to the value for which the square root argument 
becomes negative, thus the energy becomes imaginary: 

/ sin 2 F I 3a? , . [2 . . 

F ~ 7 r-ar+... => y 1 - (— + - A/ 1 - ~Y ^ l°l “ V 3 = °- 8164965809 --- (6-17) 

For a between the topological solution value a ~ —0.809466 and the minimum value 
a = —0.8164965809... F is negative at infinity, but the numerics are unstable, so it is not clear 
when it is actually reached. As an example, for F(0.01) = n — 0.008096, F'(O.Ol) = —0.8096, 
F(250) ~ —0.5 and falling (see fig 03), and for a = —0.8164965 we get F(300) flattening out 
in the neighbourhood of —25 (see fig©. However, from the equation of motion, we see that 
the allowed values at infinity are only F( oo) = kir/2, with k an integer. Of course, for odd 
k, we don’t have a solution with topological charge. 

For a larger than the topological value, F( oo) is positive, and the same remarks apply. 
For example, for a = —0.8093, F(300) ~ 0.8 and increasing (see fig®)- 

Note that from the point of view of the equations of motion there is nothing wrong with 
|a| > 0.8164965809..., just that the energy (or action) of the solution is imaginary. Then 
numerically, we can check that we obtain horizons: F becomes more and more negative, 
until at a finite point F'{tq) = —oo. For example, for a = —0.905 this happens at about 
r = 0.449, for a = —0.816502 it happens for at about r = 18.84, etc. (see figJHI). 

In conclusion, we see that Pavlovskii’s action admits topological solutions. Most solutions 
are not topological, thus unstable, and there are no solutions with horizons that have both 
real energy and good behaviour at infinity. 

One can ask what happens when we change the sign of (3 2 in the action (16.11) . We find 
that for the solution with the same asymptotics at infinity as Pavlovskii’s topological solution 
(the first two terms in the asymptotics at infinity do not depend on (3 2 ), i.e. with b ~ 7.05 
at r = 10, specifically F(10) = 0.07, F'(10) = —0.0141 (see fig®), we get a horizon at 
r = 1.8310126459463 (computer goes out of memory). Reversely, for the same asymptotics 
at zero as Pavlovskii’s solution, i.e. with a = —0.809466, for our action, F goes down to 
about 1.36 and then goes back up to a maximum slightly above 1.6, asymptoting to F( oo) = 
7t/2 = 1.5707963268 (see fig llOl) . The same qualitative behaviour (down to a minimum below 
7t/2, up to a maximum above 7t/2 and asymptoting to 7t/2) is valid for a whole range of 
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Figure 3: Topological solution 



Figure 4: Topological solution, detail 



Figure 5: Slightly below topological solution 
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Figure 6: Limit solution 



Figure 7: Slightly above topological solution 



Figure 8: Solution with imaginary energy 














values of F'(0) = a. We can easily check that F' — 0 implies F"(r 0 ) = sin{2F(ro))/T q, thus 
we have minima below F — 7r/2 and maxima above it. 

We have tried a large range of parameters a (for asymptotics at r=0) and b (for asymp¬ 
totics at r = oo), but always the same result: If the asymptotics at zero are good, then the 
solution tends to 7 t/ 2 instead of 0 at infinity, thus being non-topological, i.e. unstable. If 
the asymptotics at infinity are good, then the solution develops a horizon at finite r (F'(r) is 
infinite, and within computer accuracy, F seems finite, thus having a square root singularity, 
as expected). A detail of the square root horizon, for F(r —> oo) ~ 0.002/r is given in figllll 
We have tried values of b (at r=10) varying between 0.02 and 7.05. A few examples: b=0.02 
has horizon at about r=0.292, b=l at about r=1.044, b=2 at about 1.290, b=4 at about 
1.574. 

Let us now check under which conditions can we have a horizon, with F(r) finite but 
F'(r) = oo at a finite r = rq. From the equations of motion one finds 


F(r ~ rq) ~ F(rq) + 


'2(sin 2 F(rq) — (3 2 r\) 


r i 


(r - ri) + ... 


(6.18) 


We see that if /3 2 > 0 (Pavlovskii’s case), the solutions with good asymptotics at infinity 
indeed cannot have horizons. For them (sin 2 F(rq) — /? 2 r 2 ) < 0 at infinity and it becomes equal 
to zero at the point where the solution becomes indeterminate. Horizons can be obtained 
if r < rq, i.e. F'(rq) = oo from below, thus for solutions that have good asymptotics at 
zero. Moreover, for this we needed solutions with imaginary energy (or action). Exactly like 
for the single DBI scalar whose solutions we described analytically in section 3, these are 
actually solutions to the action with (3 2 < 0, analytically continued through their horizons. 

Indeed, we see now explicitly that in the (3 2 < 0 case, any solution with good asymptotics 
at infinity will have a horizon at some 7q, exactly like in the case of the catenoids studied in 
section 3. As mentioned, if we analytically continue these solutions past their horizons we 
obtain solutions for the (5 2 > 0 case with imaginary energy. 

We will see what happens to the temperature of the horizon with respect to the calcu¬ 
lation in the previous section. We take a more general case of a Lagrangean of the type 

\J /($) + (W) 2 — < i >2 such that near the horizon rq we have <f )/ ~ ■\//( < f>)r 1 /(2(r — rq)), as 
is the case here. Then we can easily calculate the fluctuation equation and find that near the 

horizon it looks the same way as in section 5, with \J 1 + (V<F) 2 replaced by \J /($) + (V<f>) 2 , 
but now we also get a “mass term” for the fluctuation, with 


2 /"(*) 
m = 


cm )) 5 


2 4 (/($) + (V<f>) 2 ) 

Then near the horizon 

2 / 2 
V /C — 


+ V /(^>) + (V$) 2 V[ 


V<f> 


2 (/($) + (V4>) 2 ) 3 /2 


$ 


/2 


rq 


/($) + <f> /2 4(r — 7q)+7q 


rm (6.i9) 


( 6 . 20 ) 


and we get the same result for the temperature as in the single scalar case, namely T cx (rq) 1 
(up to an infinite factor of the energy density), just that rq is now different. 
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Figure 9: Solution with horizon (/ 3 2 < 0); asymptotics at infinity of topological solution 



Figure 10: The (3 2 < 0 case; asymptotics at zero of topological solution 

In conclusion, for /3 2 > 0 (Pavlovskii’s case) we have a topological solution, but no solu¬ 
tions with horizons, whereas for /3 2 < 0 all solutions with good asymptotics at infinity have 
horizons, and there are no topological solutions. But we know from physical considerations 
that we would like an effective field theory for QCD that gives both. 

7 Proposed effective DBI action for QCD and “Skyr- 
Blon” solution for nucleons 

QCD effective action 

In fact, the action for a probe U(l) D-brane in the gravity dual geometry, representing the 
dynamics of the IR cut-off of the geometry, has both a scalar and an electric potential, thus 
having both signs in the square root. As we saw in section 3, this action admits both Blon 
solutions, that are localized lumps of electric field analogous to the Skyrme-like soliton in 
the previous section, and catenoid solutions, that have thermal horizons, whose temperature 
we have calculated in section 5. The Blons have arbitrary charges (given by their delta 
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Figure 11: Detail of horizon 


function source), but in the quantum theory we expect those to be quantized. Their analog 
for Pavlovskii’s SU(2) theory, the topological soliton, has already quantized charges, due to 
topology. 

The fact that the IR behaviour of the gravity dual can be encoded in a D-brane probe 
with a U(l) gauge field on it is not a priori obvious. However, this procedure has been used 
before to calculate meson masses. 

The usual procedure for having topological solitons of Skyrme type representing the 
nucleons is to write a nonlinear extension of the sigma model action for the SU(2) pions. 
But then as we saw there seems to be impossible to have solutions with horizons as well, as 
we know we should, from the gravity dual picture. An action similar to the one of the U(l) 
D-brane probe would give the correct physics, but why? We have neglected an important 
fact, namely that the action we are interested in has to be valid at energies beyond A qcd- 

We know that in the analysis of QCD confinement one doesn’t hope that the pions are the 
only effective degrees of freedom. They should be the only relevant degrees of freedom only 
at large enough distances (comparable with, or larger than 1 /m^). ’t Hooft ]2S] proposed an 
idea that is now the standard way to understand confinement, the “dual superconductor”. 
Out of the SU(3) QCD action one, one uses a partial gauge condition that keeps only the 
abelian projection, to the maximal abelian subgroup, U( 1) x U( 1). The normal vacuum of 
QCD is then a “dual superconductor” background: after a duality transformation on the 
abelian subgroup U( 1) x U( 1) one is in a superconductor background, where the dual gauge 
field (monopole field in QCD) is in the broken (massive) phase, due to interaction with a 
dual Higgs field, corresponding to a monopolc-antimonopole pair condensate of QCD. This 
dual Higgs field is of course also massive in the confining vacuum. 

Then a quark-antiquark pair of QCD (forming a meson) will serve as endpoints for tubes 
of confined electric flux, as in a dual type II superconductor, giving a linear potential. In the 
dual abelian theory, monopole-antimonopole pairs (solitons of that U( 1) x U( 1) theory) will 
have between them magnetic field flux lines confined to tubes. In this dual theory, the gauge 
field will be nonzero inside the flux tube and zero outside (where it is massive), and the 
Higgs held will be zero at the core of the flux tube and go to the symmetry-breaking value 
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on the outside. Nucleons will have three quarks, thus are harder to explain, and one has 
to invoke some spherically-symmetric version of the meson mechanism. Also, in this picture 
the quarks are just sources, but if they are fundamental, as in real QCD, they introduce 
extra light degrees of freedom, the SU(2) pions (in fact, the lightest). Also, the confinement 
picture changes, as quark-antiquark pairs can break off an electric flux line. 

Seiberg and Witten m provided the first example of such a mechanism, for the SU(2) 
J\f = 2 SYM theory broken to the confining J\f = 1 SYM by a soft-breaking susy mass 
term mTr$ 2 . The low energy effective action for the SU(2) J\f = 2 SYM theory is an 
abelian (U(l)) J\f = 2 susy theory. At the point on the moduli space where monopoles 
become massless, one has to make a duality transformation and go to a description in terms 
of monopoles. When one adds the susy-breaking mass term, one is driven to the (former 
moduli space) point where monopoles are massless, and one finds that the dual Higgs held 
(monopole-antimonopole held) condenses, making the dual photon massive. This is the 
picture expected from the ’t Hooft idea, just that it is valid only at low energies (on the 
moduli space). Seiberg and Witten calculate exactly only the low energy theory, given by 
f d 4 xd 4 0J r ( ] &), but higher order corrections (coming from terms like f d 4 xd 4 6d 4 67i(^’ 
for instance) are not calculated. 

So for the effective theory of QCD in the confining vacuum we expect to have a dual 
abelian (probably U( 1) x U( 1)) theory for a gauge held interacting with a Higgs that gives 
the photons mass, and also the SU (2) pions (almost Goldstone bosons for the chiral symmetry 
introduced with the light fundamental quarks). The low energy U( 1) x U( 1) gauge theory 
is probably complicated (as seen in the Seiberg-Witten example), and the low energy SU(2) 
pion action is the usual nonlinear sigma model action 

S = I d 4 x^[TrL^L fl + mlTr(U + U + - 2)] (7.1) 

where as before L M = U~ 1 d[ J U. At high energies we don’t have much information, except 
from the gravity dual theory. The action for a single scalar pion was of DB1 type, as it 
corresponded to the position of a brane (IR cut-off) in an extra dimension. We know that 
the QCD pions are an SU(2) triplet, probably corresponding in the dual gravity theory to 
fluctuations of the IR cut-off (physical brane) in an SU(2)-valued set of directions, like on 
an S 3 . Thus we want to have an SU(2) generalization of the pion DBI action that reduces 
to ED at low energies. 

For the dual gauge degrees of freedom, we will take a single U(l) (probably the correct 
theory involves both U(l)’s, but we take the minimal assumption). We will also assume that 
for the problems at hand, the dual Higgs will only give mass to the dual photon. If the 
dual superconductor is correct in all details, the Higgs should change from zero inside the 
nucleons to the symmetry breaking value outside. However, from the gravity dual point of 
view an action for the Higgs would imply making assumptions about the actual model for 
gravity dual of QCD. Since we don’t have such a model, we don’t have any idea how such 
an action would look like. The minimal assumption then is to hope that the dynamics of 
the Higgs is irrelevant for the nucleons, and either it is approximately constant, or at the 
surface of the nucleon, where the Higgs would vary, the pion held is more important than 


32 


the Higgs. In that case, we could just replace the Higgs with a mass term for the U(l) gauge 
field. But what would be the high energy behaviour of the U(l) gauge action? From the 
dual picture, if we can think of the U(l) as the gauge field on a probe D-brane, the action 
should also be of DBI type. It is also hard to see what else could it be, as the DBI action is 
the unique nonlinear correction of the Maxwell theory that is both causal and has only one 
characteristic surface (generically, there are two) [28] . 

However, that can’t be the gauge action for arbitrary energies either, we expect that 
above some energy the simplest description of the gauge fields should be in terms of the 
original gluons. Thus this action should be valid only for U(l) energies smaller than some 
scale M, above which we cannot say anything. The simplest assumption that we will use 
afterwards, is that effectively |V0| < M, i.e. after it reaches M it stays fixed at that value, 
but in reality it should be a more complicated condition. The gravity dual picture suggests 
something similar, as we will discuss in more detail later. Indeed, we have seen that at Er 
in the gravity dual the created black holes start feeling the curvature of space. Thus at this 
energy gravity (dual to the gauge fields) starts feeling new terms in the action, in agreement 
with having a relation between M and E R . 

Thus we propose the action (at zero magnetic field) 


S = 


flMl 


d A x 


p 2 r 2 


{EiLiY 


LI 


m 


+ 772 + 772 V + U*-2) + - 1] (7.2) 


M'tM'i ' Ml Ml 


where E = V0 and Ma is the gauge field mass, and E max = M. 

The most conservative ansatz for the mass scales is M\ = M 2 = Ma = A-qcd, but they 
need to be only of the same order of magnitude from physical considerations. Most of the 
time we will assume this ansatz to be true however. 

A priori there should be some Chern-Simons term in the action as well, that reduces to 
the one in the original Skyrme model [IT2] . and since also the D-brane action in the gravity 
dual will have in general such a term. But we know that it will not contribute for spherically 
symmetric solutions (which we are interested in), so we will ignore it in the following. 

D-brane action toy model 

To understand the issues better, let us look at the simpler model of the single scalar pion, 
described by the U(l) D-brane action. This is the case analyzed in the gravity dual, with 
the D-brane scalar being the position of the IR cut-off (brane) in AdS space. Moreover, we 
analyzed its solutions in section 3, so now we can interpret them. In section 3 we looked 
at the massless action, which had as solutions the Blon, the catenoid, the BPS Blon and 
general solutions interpolating between these three. 

In the gravity dual, the cleanest case is the case when the scalar is also absent (pure 
gravity, the IR cut-off is fixed, i.e. the radion has a very large mass). That would correspond 
to putting the DBI scalar to zero. Then the Blon should stand for a glueball-type solution 
(excitation of the theory). The constant C in the Blon solution can be taken as electric charge 
(we identify it as such by the behaviour at infinity), thus is quantized, as it should. But if 
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the gauge field is massive, as we have argued it has to be (being in a dual superconductor 
background), we will not measure an electric charge at infinity. However we also know that 
by colliding two glueballs we create black holes in the gravity dual, and eventually these black 
holes live on the IR cut-off, being effectively 4 dimensional. Thus eventually (in the Froissart 
regime) in the 4 dimensional collision of two field theory glueballs at very high energy we also 
expect to create a fireball type solution with a horizon. But we saw that the pure electric 
action (no scalar) has no solutions with horizons. It must follow that the description of the 
pure Yang Mills case at large energies cannot be in terms of the pure U(l) DBI action, but 
must include a true scalar (for the description of a scalar glueball field), which scalar can 
generate solutions with horizons. This scalar should be relevant above some energy scale Ep 
(the unknown Froissart energy scale for glueballs). It could be for instance that the Higgs 
field giving mass to the U(l) (which we ignored in this analysis) has also a DBI-likc kinetic 
term with energy scale Ep. But there is one more modification that one expects to the U(l) 
DBI action: at energies above the soft Pomeron scale Ep, in the gravity dual the created 
black holes start feeling the curvature of the space, thus the effective action for U(l) should 
get new terms. They will modify the value of (j) only very slightly, but enough to ensure that 
for instance the energy density at the center of the Blon doesn’t diverge anymore, but stays 
finite (the energy density is proportional to 1/yR — (j)' 2 ). 

In the case that we have a nonzero single scalar pion (in the gravity dual the radion 
has a small mass), the Blon would again correspond to a glueball solution. Now a Blon-like 
general solution with C > C will be also a type of dressed glueball, and by colliding two such 
solutions (at least for two solutions with opposite charge C so that the result of the collision 
is neutral) one should create a scalar pion fireball, with a horizon. Such a solution is the 
catenoid, and catenoid-like solutions with C > C if the total charge of the colliding objects is 
nonzero. But now there should be not only glueballs in the theory, but also hadrons, which 
have a scalar pion profile and quantized scalar charge (toy version of the baryon charge). We 
need then C > C to avoid horizons, and if both have unit values, C — C — 1, which would 
imply the BPS Blon is such a hadron. However, we don’t want for hadrons solutions with 
diverging scalar at zero, so we can either consider C > C and find a solution with a given 
finite value of the scalar at the origin X (0) = X 0 (which specifies C if C=l), or consider 
that the action has higher order corrections that limit the value of the derivative JY(0), such 
that the BPS Blon becomes also finite. In this toy model, the baryon charge is not well 
defined, but we will see that in the real case we will have a combination of the two proposed 
solutions, namely we will have a given value at zero A"(0) and the scalar derivative A'(0) 
will be limited in value. 

Solutions and SkyrBIon as nucleon 

We will first neglect the mass terms (m 2 and m\ terms) in (17. 2 j) . Note first that for 
E = 0 we get the Pavlovskii action with /3 2 < 0, case we analyzed in section 6, and for U=1 
we get the original action of Born and Infeld. Since the fields appear quadratically in the 
action, the above truncations are consistent, and we can embed any solutions we found in 
those cases in the action (H with = Ma = 0. 

We can easily check that (njLj) 2 = —F' 2 , thus for a spherically symmetric solution, for 
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which E % = —di<p(r) ~ ?+, ( E{Li ) 2 = —(j) ,2 F' 2 and the action on the solution becomes 


S = — I dt r 2 dr[\l (1 — 0 /2 )( 1 + F' 2 + 2 ^-—^-) + 0 /2 F /2 — 1 ] 


dt r 2 dr[\l (1 - 0' 2 )( 1 + 2^-^) + F /2 - 1] 


with equations of motion 

r 2 0' + 2 (j)' sin 2 F , 

^ R J 


= 0 ; 


r 2 F' (1 — 0 /2 ) sin 2F 


R 


R 


where R is the square root in the action. 

The first equation can be easily integrated with a constant, giving 


C 


R = 


C 2 + r\ 1 + 2^) 

r 4( 1+2 


'1 + 


F' 2 


1 + 2 sF^F 


1 + F' 2 + 2 


sin 2 F 


(7.3) 


(7.4) 


(7.5) 


C 2 + r 4 (l + 2+£) 

Then the second equation is (16.711 (for /3 2 = —1/2) with a modified right hand side, i.e 


(r 2 + 2 sin 2 F)F" + (2 rF' — sin 2 F) 

3 1 

+2(rF /3 - F' 2 sin 2F + -F' sin 2 F -- sin 2 F sin 2 F) = 

/y* ry* £ 

1 + F' 2 + 2 sin 2 F/r 2 C 2 

l + 2sin 2 F/r 2 C 2 + r 4 (l + 2 sin 2 F/r 2 ) 

sin 2 F 

x[2rF' — sin2F + 2sin 2 F( --—)] (7.6) 

ry > ry* £ 

We can now find the behaviour of solutions at r = 0 and r = oo. Let’s first assume that 
F( 0) = Ntt is a good vacuum at r = 0 and expand around it. Defining </o = 0(0) = 0'(O) 
and F 0 = F( 0) = F'{ 0), we get from (17.511 that 



i+3 n 

(1 + 2 FS) 


(7.7) 


Note that this is the condition that R, the square root in the action, is zero at r=0, the same 
condition we obtained for the U(l) D-brane action in section 3, in particular for the Blon 
solution. Thus as for the DBI action studied in section 3, 0 O = 1, F 0 = 0 is a solution, and 
is actually the same solution, the Blon, since for F=0 we have the same action, and this is 
a consistent truncation of the original action. 
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However, when we try to satisfy fLU at nonzero F, we see that it is impossible. A Taylor 
expansion doesn’t work, as the 0(r) terms on the left hand side of (17.61) cancel, whereas on 
the right hand side we get 


1 + 3F 0 2 

1 + 2Fq 


[—2 rF*} 


(7.8) 


so that would imply F 0 = 0. But then we would get 


(r 2 F")(0) = -2(rF' 3 )(0) (7.9) 

which doesn’t have any solution for any nonzero coefficient in the Taylor expansion. A 
singular behaviour of the type F — Nn = ar a , with 0 < a < 1 doesn’t work either, as it will 
give the equation 

a 3 r 3a ~ 2 a(a - 1 + iV3)(a - 1 - iV3 ) = 0 (7.10) 

with no real solution between 0 and 1. The other possible singular behaviour, F' — a In r, F — 
Ntt ~ ar (In r — 1) doesn’t give a solution either, the left hand side being of the order rln 2 r 
and the right hand side of the order r In 3 r. 

So if we would only consider m we would obtain that, as in section 3, 0o — 1 => Fq — 0, 
but now moreover 3/2 > 0 q > 1, and 0o increases monotonically with Fq. Based on the 
toy model of section 3, we would expect to find solutions corresponding to all these values, 
however as we saw, m implies that in fact that all solutions must blow up before reaching 
r=0, as in the single scalar case! (Even in the pure pion case, at 0 = 0, we had solutions 
that went to F — Nn , though they had the wrong asymptotics at infinity, but now even 
those are excluded). 

For the behaviour at r = oo, we find as before, that the noninteracting asymptotics is 
not modified: Putting (and assuming F( oo) = 0) 


</> 




a 

r^n 7 



(7.11) 


we find n = 2, m = 3 as for the case where the scalars don’t interact with each other. 

For 0 = 0 we have as we mentioned Pavlovskii’s action for /3 2 < 0, for which we saw that 
any solution that has good asymptotics at infinity will have a horizon. 

For a general solution (at nonzero 0) with a good behaviour at infinity for F, let’s under¬ 
stand the formation of the horizon. If the horizon would be only in F, i.e. F' = oo, F finite 
at r finite, and 0' finite, we would get 


F ~ Ft 



(rf + 2 sin 2 F(r 1 ))(l — 0 /2 ) 


r i 


n) 


(7.12) 


which is equal at 0 7 = 0 to the condition we had for the Pavlovskii case. However, in fact 
we can see from (H3D that if we have F finite, F' infinite at r finite, we also have 0' infinite. 
Then we get 


F ~ Fr 



r\ + 2 sin 2 F(r\)) 


n(l - a) 


n) 


(7.13) 
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Figure 12: Horizon detail 


where 

C 2 1 + sin 2 F\rq)/r 2 

C 2 + rf(l + 2 sin 2 F(ri)/rf) 1 + 2 sin 2 F{ri)/r\ 


(7.14) 


thus we always have a solution, as in the Pavlovskii case. 

We also deduce that there can be no horizons of this type at r=0, where 0' < \l~3~j2 
(with F finite and F’ infinite, since that would need 4>o <1, which is excluded), as we have 
already argued. 

Finally, we have done an extensive numerical solution search for the equation (ED con¬ 
firming that there are no solutions that go to F = Nix at zero, and the solutions with good 
asymptotics at infinity have horizons at finite r. As the parameter b in the asymptotics 
at infinity of F in (I7TT1) is decreased, the position of the horizon rq decreases, and F in¬ 
creases until a maximum around F ~ 1.478 for b ~ 0.29 and rq ~ 0.172 (all this choosing 
C=l; see the horizon detail of this solution in fief II2D . An example of horizon solution with 
b = 0.6, rq ~ 0.48 and F m ~ 1.31 is given in figTHl and an example with b = 0.02, rq ~ 0.001 
and F m ~ 0.28 us given in fig HTl both of them having C=l. As one increases C, the 
maximum of F increases as well: for instance at C=3, for b ~ 1.8 we have F ~ 2.18 at 
7q ~ 0.36. 

The conclusion is that this action (O does not have topological solutions, contrary to 
the intuition that we gained in the single pion case, the toy model D-brane action. However, 
we have to remember that we said the action is only valid up to \E max \ = M. At F=0 
there is no problem, as then \E\ is bounded by M 2 , but in the presence of diverging F', \E\ 
would diverge as well. 

As we mentioned, the simplest assumption is to say \E\ becomes constant and equal to 
M (iu reality it could increase at a slower rate and saturate at a higher energy value). But 
by the relation between F' and <f>', F' will be forced to saturate as well. Indeed, then the first 
equation in Q can’t be easily integrated to give m anymore, but the modified second 
equation can be seen to still imply that F’ is infinite only if 0' is infinite. We observe from 
m, valid before saturation, that the closer we are to r=0 when 0' saturates, the larger the 
saturation value of F'. 
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Figure 13: Solution with larger b than the maximum 




Figure 14: Solution with smaller b than the maximum 
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Finally, this means that we will nevertheless have a topological solution, reaching F = n 
because instead of the horizon in F at finite r of the action OD, we can continue with 
> M down to r=0. Numerically, this will work if M ~ 10 (in units of Mi = M 2 ), 
since as we saw, for the maximum, 7 r — F max ~ 1.7 and rq ~ 0.17 (and moreover, for larger 
r 1 , F max decreases only slightly). Without a complete high energy action, we can’t derive 
the exact solution unfortunately. 

We will call this solution the skyrBIon and it will represent a nucleon. 

Before we analyze it further, let us look at a point that may be of concern. Since we 
avoided the horizon in F at small r, why do we still have solutions with horizons at all (as 
we claimed, the solutions with horizons correspond to fireballs created in collisions)? The 
answer is that we have treated until now the massless case. But of course, we have nonzero 
masses for the pion and the gauge held. The mass of the gauge held should be around the 
hrst glueball mass, thus around M|, M 2 . The mass term will not affect the solutions at 
r < M a rMi ~ 1 (in our units), but at larger distances, 0 will decreases exponentially to a 
negligible value. Thus if F' diverges at r > 1 it will not yet be coupled to the (very small) 0 
held and thus F' will not be bounded. Only at r < 1 as for the skyrBIon solution this can 
happen. 

Interpretation 

Since the proposed skyrBIon solution will represent a nucleon, it has to have the right 
properties. The topological charge it carries, B, will be interpreted as baryon number, as 
usual. It can also be positive or negative, corresponding to particles and antiparticles, as 
usual. It is not clear from what we derived whether there are solutions with higher baryon 
number, that will depend crucially on the high energy modification of the action. What 
seems certain is that there aren’t any solutions with arbitrarily high baryon number, there 
will be a maximum B. That is so if the slope F’ cannot increase indefinitely, since the point 
r\ at which it starts saturating cannot be larger than 1/A qcd- That is in accordance with 
the real world, where stable particles (spherically symmetric configurations) of arbitrarily 
high B don’t exist, and solutions of B > 1 if present could represent possible metastable 
particles. This situation is to be contrasted with, for instance, the case of SU(2) DBI pions 
of “wrong” sign studied by Pavlovskii, where an arbitrarily high B was possible. We should 
mention that by contrast, nuclei (which have higher B), will not be spherically symmetric 
configurations, as is true in the real world. We will discuss them shortly. 

The skyrBIons have also a (quantized) U(l) electric charge. Due to the mass of the 
gauge held, it is not a charge detectable at infinity, but we will have skyrBIons of negative 
U(l) charge with the same (0, F) profiles (the equations of motion are invariant under 0 
reflection), but different U(l) interactions. But that just accounts for the fact that the 
proton is different than the neutron, even once you take the electromagnetic interactions 
out! I 11 other words, p-p vs. p-n scattering differs more than by the purely electromagnetic 
exchange diagrams. So the massive U(l) charge +1 and -1 particles with B=+l are the p 
and the n. 

We also have the Blons in the theory, purely U(l) solutions with an arbitrary quantized 
charge, but these are easily seen to be undestood as scalar glueball excitations, since the 
massive U(l) should be dual to the SU(3) gauge fields. 
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Nuclei as Bionic crystals 

Finally, how do we undestand the nuclei in this picture? It was understood for a long 
time that in Skyrme-like models, the nuclei should be bound states of nucleons forming sort 
of a crystalline structure and not by a higher B spherically symmetric soliton. This is indeed 
what we obtain, as we shall see, with the added bonus that it seems impossible to even have 
arbitrarily high B solitons anyway. 

First, notice that Blon solutions of the massless BI action (U(l) gauge held) can form 
crystals (221 iSf of the NaCl type, i.e. with positive and negative U(l) charges alternating 
in the crystal. This is maybe not so surprising, but one can write down an exact solution 
for this. If the Blons in our theory represent glueballs as argued, it would be interesting to 
see what this fact could imply experimentally for glueballs. Of course, as we said, the gauge 
held is actually massive, so the Blon crystal would be modified in our theory anyway. 

More importantly, let’s see what will happen for the skyrBIon solution. They come 
with massive U(l) charges +1 and -1, that we argued should correspond to the proton and 
neutron. At large distances (r > Ma) only the pion interaction is relevant, and it is an 
attractive interaction, being mediated by perturbative scalars. That is, a two skyrBIon 
solution at large distance will have a smaller energy than two individual skyrBIons. 

As we argued at the end of section 3, at small distances, as long as the energy is finite, 
we can hgure out if the potential is attractive or repulsive from the behaviour as a function 
of charge. The energy of the static spherically symmetric system in the case of zero mass 
and no higher order corrections is 


H 



1 + F' 2 + 2 sin 2 F/r 2 
R 


i] 


(7.15) 


We see that only the sin 2 F/r 2 terms stop the action from scaling as charge to the 3/2 power 
as in section 3, and if the radius where F’ and 0' reach their maximum is not too small, this 
term will not affect too much the scaling. We also see that before even reaching radii smaller 
than 1/Ma ~ 1/A qcd and the U(l) held becoming relevant, the pion scalar held becomes 
nonlinear, and given the analysis in section 3 could already become repulsive. We can thus 
say that the interaction will certainly turn repulsive at some distance inside 1/A qcd, maybe 
even before. 

Finally, we can say that once we reach the zone of saturation of the derivative of 0 
(at r ~ ri) the potential will be clearly repulsive. Indeed, it is probably safe to assume 
that the energy density will be approximately constant (p 0 ) inside it (for r < r\). If F’ is 
approximately constant inside r \, then r\ is proportional to the topological charge B=N (at 
r=0, F = Nir), thus the energy in the core will increase like p^rf ~ N 3 , thus two unit charge 
baryons will be strongly repulsed. 

Given this picture it is clear that there will be a minimum rc for the potential between 
two nucleons that is at least for radii inside 1/A qcd < 0.2/m if not outside it (larger). 
Moreover, we talked about the interaction of two identical nucleons, but there will clearly be 
a diffence between the U(l) like charge and opposite charge interactions at least at distances 
of the order of r\. 

All of the above features are in accordance with the real world, where a meson (pion) 
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interaction potential works down to about 0.4fm after which the interaction potential rises 
dramatically. 

The structure of nuclei is known to be well described by the Bethe-Weizsacker formula 
for the binding energy of the nucleus 

B(A, Z ) = a v A - a s A 2 ' 3 - a c + a P 5A ~ 3 / 4 - a sy J N (7.16) 

where 5 = +1 for (N,Z) even-even, -1 for odd-odd and 0 for even-odd or odd-even. The 
first three terms (volume term, surface tension term and Coulomb interaction) are known to 
be deduced from a liquid drop model, and the a P term is quantum in nature (the nucleons 
are fermions of spin 1/2). For the liquid drop model a nucleon-nucleon potential with a 
minimum is probably enough, and that was qualitatively derived already. The last term in 
the binding energy can also be understood qualitatively in our model. As stated before, in 
the U(l) D-brane toy model, the Blons form NaCl-like crystals, with alternating +1 and -1 
charges. Thus the minimum energy in this case is given by equal number of +1 charges and 
-1 charges. The same should be true if we add some (constant, positive) scalar charge to all 
Blons, such that we keep C > C. 

It is then very likely that for skyrBIons also the minimum energy configuration is given 
by equal numbers of U(l) charges +1 and -1, that we equated with n and p, thus for N=Z. 
Thus we expect indeed a term that drops with (N — Z) 2 in the binding energy (the fact that 
it is (N — Z) 2 and not some other power is just because a perturbative U(l) interaction is 
proportional to the charge squared). 

8 High energy scattering models 

In this section we will see what we can derive about the high energy scattering of hadrons 
that are modelled by boosted field theory solutions. We will look both at the toy model 
solitons and at the skyrBIons. 

We will start with a few general observations. If we take two boosted hadrons, represented 
by boosted solitons, and collide them, we can follow Heisenberg’s argument and find again 
the cross section: 

If we have a coefficient of inelasticity (energy loss is £ — ay/s) 

( 8 .i) 

then, provided that the action is of DBI type, Heisenberg showed that b max is obtained when 
£ =< E 0 >~ m 7r (for a linear action he gets < E 0 >~ 7 ~ /s, thus gets a b max independent 
of s, so it’s not relevant for the high energy behaviour of real hadrons). 

Then Heisenberg suggests that a should be proportional to the pion wavefunction overlap, 
and the wavefunction is proportional to e~ mr , thus a ~ e~ mb : 

j d 3 rM'r)Mr) ~ J d 3 re~ mr e~ m{b - r) ~ Ve~ mb (8.2) 
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after which one obtains a ~ ln 2 s. In our case though, if a ~ 1/fP, then we find cr ~ s 1//p . 
But now if a is due to a wavefunction, pion or otherwise, and it satisfies i(j(r) ~ 1 jr n , then 


a ~ / d 3 r'0i(r)'0 2 (^) ~ / d 3 r- 


ry* I o I ^ 


6 |' 


})2n —3 


(8.3) 


thus p = 2n — 3, contrasted with the case when ex ~ ^(b) which gives a ~ l/r n , thus p = n. 
Note that in Heisenberg’s case there is no difference in between the two hypotheses about 
a, nor if we also take derivatives on the wavefunctions, in all cases we get a ~ e~ mb . Now 
however, we get different results. Moreover, there is no good reason why the wavefunction 
overlap behaviour with s would dominate over the dynamics of the theory (in the Froissart 
regime, the overlap is exponential, thus dominates). 

In the gravity dual case, we had a A-S shockwave perturbation in flat space or AdSd+i x 
Xj, with 

|V*M| ~ £ (8.4) 

and the horizon was situated where |VT(r)| = 1, giving b max ~ s 1 ^ 2p \ and since a ~ 
nb 2 max ~ s 1/p , it is the same kind of analysis as we have now, if we identify |V\l/(r)| with 
S(r) = y/scn(r) (it would be also with y/s^){r) if a VO- 

We should note that the quantity of interest is y/sa(r), which increases with y r s. The 
pion held is a scalar, thus as we noted - 0 (r) (and thus a{r ) when it is only due to the pion 
wavefunctions) stays constant as one boosts. 

From the collision of high energy boosted solitons (Blons or skyrBIons) it is hard to 
derive the behaviour of ct(r), but we can try to understand at least the energy regimes. All 
energy regimes of interest obey yfs M ~ A qcd and then the soliton becomes a shockwave 
(corresponding in the dual to yfs > Mp , when particles become gravitational shockwaves). 

Now we will adress specific cases. First, when there are no pions at all (in the gravity 
dual we only have gravity, the IR cut-off is non-dynamical), as we argued in the last section, 
we should have the massive U(l) held but in order to generate horizons we should have also 
a scalar held becoming relevant at the Froissart scale Ep. We will assume for generality that 
the mass of the gauge held, Ma, (as well as the mass of the new scalar held, M' A , that should 
be of the same order) is smaller than the DBI scale of the gauge held, Mi. Then we have a 
region 1/M A > r > 1/M\ where we can still neglect the masses. 

We have seen in section 5 that if we can ignore the masses of the helds, the DBI action 
has a solution with a horizon and a temperature that goes like T ~ M -1 / 3 . In the collision 
of two glueballs we will create such a scalar held (Higgs?) fireball, thus there should be some 
evidence of thermalization in the final state (after the decay of the fireball). However, this 
object will not be dual to a black hole (but rather to black hole creation integrated over the 
5th dimension), so its analysis will be more complicated. It is hard to estimate the behaviour 
of a(r ) and thus of er(s) at this point. It should depend on the massive gauge held prohle 
and on the scalar held (Higgs?) prohle. As one continues increasing the energy, eventually 
the core of the Blon will become relevant. As we argued, when the scattering energy reaches 
the “soft Pomeron scale” Ep, in the gravity dual the produced black holes start feeling the 
curvature of the space, and correspondingly the U(l) DBI action should have high energy 
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corrections that keep the energy density from diverging at the core of the Blon. As one 
boosts the Blons above Er, part of the interaction is due to the overlap of the Blon core 
with the tail of the second Blon, 0 ~ Cifr. Thus the high energy corrections to the Blon 
become important and the scattering cross section behaviour should be modified too. When 
r 0 ~ 1 /M' a , thus when y/s ~ M F /M A 3 , one will still create a fireball, but the temperature 
will be different. The cross section behaviour however still depends on the dynamics of the 
theory, at least until an unknown energy scale Ef • Note that Ep > Er > M F /M A 3 , since 
from the gravity dual we know that Er = M f / Aq CD = N 2 Aqcd > Mp/ Aq CD an d Ef > Er. 
Thus the temperature of the fireball will be modified before the onset of the modification 
of the cr(s) behaviour (which should occur at Er > M F /M' A ). Finally at Ep we will be in 
the Froissart regime. Indeed, a perturbative calculation for the horizon of the formed Blon 
at tq 1 jm-K (how much larger depends on what Ep is) gives 0 ~ Ce~ MAr /r ~ 1 , thus 
b ~ l/M A \n(CM A ) and C oc M a (a=number). 

If we have a single scalar pion, i.e. the U(l) D-brane toy model (in the gravity dual the 
IR cut-off is dynamical now), we will also have a pion wavefunction. The fireball that will be 
created in the scattering will now be due to the pion held, which dominates over the second 
scalar, having larger wavefunction. The behaviour of a{r) and thus of cr(s) until r$ 1 /M' a 
will still be dominated by glueball physics, but now the onset of Froissart behaviour is 
sometime after r 0 ~ 1 /m^, thus sometime after -y/s ~ Mp/m 3 , since the fireball is made up 
of the pion scalar held. In between ro ~ \/M' A and r o ~ 1 /m^ the pion held wavefunction 
predominates, but we still expect the dynamics of the theory (QCD) to dominate over 
the wavefunction overlap. Now we have the same onset of the soft Pomeron behaviour at 
Er = Mp/A 7 Q CD = N 2 Aqcd > Mp/Aqcd (thus when r 0 is already > 1/M A ), but the 
Froissart behaviour will onset even later than before, at E' F . After E' F (when r 0 1 /m n ) 
we have Froissart behaviur in terms of m T instead of M A . 

Finally, if we have SU(2) pions, and the nucleons are modelled by the conjectured skyr- 
Blon solutions, the overall picture is the same as for the single pion. The only difference 
is that now the skyrBIons will carry a topological charge (baryon number), and again as 
before the energy density at the core of the solution will not diverge, due to higher order 
modifications to the action. 

Another important issue to discuss is the issue of transparency of the scalar “fireball”. 
We have seen that in all cases (above Mp, even before the Froissart bound) we will produce 
metastable scalar held solutions with thermal horizons. We have also argued in section 5 
that the horizon of the fireball acts like the horizon of a black hole with respect to scalar 
excitations, i.e. it will infinitely delay the exchange of information with the outside. This is 
exactly true if the temperature of the fireball is nonzero and finite, and if the temperature 
is infinite, is is true only for high energy modes (the phase and group velocities go like 
l/Vk). In )B] we have argued that the “jet quenching” observed at RHIC, which we argued 
should be in the Froissart regime, is nothing but the information paradox of black holes, 
i.e. information (quantum particles) coming in, and thermal radiation coming out. But 
we know that “jet quenching” is not absolute, there is information coming out, even in the 
RHIC (Froissart) regime, and certainly before that, so how do we reconcile this with the 
observation made that fireballs have information absorbing horizons? 
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In the gravity dual, it is easy to understand this. One integrates over the fifth coordinate, 
i.e. the black hole produced in the bulk has quantum fluctuations along the fifth dimension. 
As one increases the energy, the black hole becomes more classical and gets stuck on the 
IR brane. Thus particles thrown at the interacting region in 4d held theory can miss the 
black hole in the fifth dimension due to quantum fluctuations, but that becomes incresingly 
rare as the energy of the produced black hole increases. But how do we understand this 
phenomenon in held theory? 

First oh, we saw that for the catenoid (classical, static, spherically symmetric solution) 
the temperature is infinite and then the velocities at the horizon go like 1 /Vk, thus the 
infinite information delay at the horizon is only approximate, and valid only for high energy 
modes. Since the RHIC experiments deals with high energy probes anyway (jets, hence the 
name “jet quenching”), this is exacly what is observed. But we argued that the temperature 
becomes finite, due either to time dependence, higher order corrections to the action or 
quantum effects, in which case the horizon should absorb all information. But while it is 
true that the fireball horizon absorbs scalar held information, it is also essentially (i.e. even 
classically) unstable: one has to continue the scalar held solution inside the horizon, and 
the most likely possibility, the “brane bubble” depicted in hgJ3 is unstable, and there seems 
to be no way to continue to a singularity behind the horizon. Another symptom of the 
classical instability is the infinite temperature of the simple catenoid, for which we said a 
possible cure would be through time dependence. The fireball also has a horizon, so it is 
clearly quantum mechanically unstable also due to the thermal radiation. But unlike the 
black hole, the absence of a singularity that would crush everything makes it unclear why the 
information would be destroyed in the first place. A scalar particle entering the metastable 
scalar fireball would linger at the horizon for a long time, but as the fireball horizon would 
dissappear it could presumably continue forward. If the time scale of the fireball existence 
is long enough, the particle could thermalize, giving rise to “jet quenching”. 

Note therefore that the scalar fireball is a much cleaner example of the information 
paradox, meaning that quantum particles collide and thermal radiation comes out, in a 
purely quantum mechanical scattering. This underscores the fact believed in string theory, 
that the black hole really doesn’t destroy information G3, and the information can be 
retrieved from the almost thermal radiation coming out. The scalar fireball production 
is also a first step towards finding a held theory formalism where the temperature is not 
introduced by hand, but appears due to the time evolution (a zero temperature system 
creates a finite temperature one through time evolution). 

9 Discussion and conclusions 

In this paper we have analyzed the nucleons (baryons) and their high energy scattering in the 
fixed t, high s regime from the effective held theory point of view, guided by the gravity dual 
description developped in me™. At rest, the nucleons are supposed to be described by 
the Skyrme picture, as topological solitons of the pion held, and at high enough energy, by 
colliding shockwaves of the pion held, according to Heisenberg’s model |8| for the saturation 
of the Froissart bound. We wanted a description that can interpolate between these two 
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cases, with the nucleons being solitons of some effective action for QCD involving the pion 
held that when boosted to high energies become the colliding shockwaves. 

Both in the Heisenberg model and the gravity dual descriptions, the DB1 D-brane action 
played a major role, so we studied its static solutions, the Blon, the catenoid, the BPS Blon 
and the solutions interpolating between them. The D-brane action was a toy model for the 
effective action for QCD, with the solutions interpolating between the Blon and the BPS 
Blon standing for the nucleons, and the solutions interpolating between the BPS Blon and 
the catenoid standing for fireballs created in the high energy collision of nucleons. Since at 
the energy scale Ep, in the gravity dual the action is changed by the fact that the created 
black holes start feeling the curvature of space, it follows that at high enough energies the 
same should be true for the U(l) effective gauge field. Its action should be modified at high 
energies, thus ensuring for instance that the energy density at the center of the Blon doesn’t 
diverge. We analyzed the boosted solutions and then the fireballs that will be created in 
their collision. The “dumb holes” of Unruh uni are an example of held theory (specifically, 
hydrodynamics equations) solutions that have thermally emitting “horizons” where v = c 
and dv/dr is finite, such that (in Unruh’s case, when p and c are nonzero and finite at the 
horizon) 


m 1 dv 

T = -- 

2^r dr 


(9.1) 


They can be in fact mapped to the black hole solutions, hence their name. We have shown 
that catenoids are also of the same type, we can also map the fluctuation equation to a black 
hole equation, and the surface where the scalar field diverges, i.e. where X' —> oo, while X is 
finite, also has v = c and dv/dr is consistent with nonzero temperature. We calculated that 
the temperature of the static spherically symmetric solution is actually infinite, but argued 
that the infinite prefactor should be regulated, either classically or quantum mechanically 
and then the temperature scales as 


T oc 
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M 1 / 3 


(9.2) 


where r 0 is the horizon position, Mp is the DBI scale, and M is the mass of the solution. 
We similarly gave a perturbative argument for the case where the scalar field has a mass 
m, that the temperature of the modified catenoid will be asymptotically (for large mass M) 
proportional to m. The horizon of the catenoid can be probed with other types of fields, 
but scalar excitations propagate as in a black hole background, thus the propagation of 
information, described by their characteristic surfaces, will be obstructed (infinitely time 
delayed for high energy modes, the phase and group velocities scale as 1/ y/k) as for a black 
hole horizon, even though it takes a photon a finite time to go to the horizon and back 
(since we are in Minkowski space). We analyzed the propagation of waves in the scalar 
background, and proved that if the background corresponds to a black hole of finite nonzero 
temperature, then light takes an infinite geodesic time to reach the black hole horizon, and 
correspondingly for the scalar, the phase and group velocities go to zero at the horizon for 
all modes, besides scaling as 1 /y/k. 

In the real world, the pions transform under a global SU(2), giving the possibility of 
having topologically stable solutions, that one identifies with the nucleons in the Skyrme 
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program. In ra it was shown that a wrong-sign SU(2) generalization of the DBI action has 
Skyrme-like solitons. We analyzed this action in detail and found that there are no solutions 
with horizons, only topological solutions. For the action with the correct sign, we have again 
analyzed in detail and found that any solution with the correct asymptotics at infinity will 
have a horizon, and there are no topological solutions. The temperature of the horizons is 
now also T oc (rq)" 1 , but the position of the horizon, rq, is modified with respect to the 
single scalar case. 

The creation and decay of the scalar fireball was found to be a clean laboratory for 
the black hole information paradox, since the same situation seems to appear: in high 
energy quantum collisions one creates thermally radiating fireballs, apparently violating the 
unitarity of quantum mechanics. But this shows that the problem of the information paradox 
is not due to quantum gravity, but rather to the lack of a unifying held theory formalism 
where finite temperature can appear in the process of purely quantum mechanical scattering 
at high energies, as we know happens in the RHIC experiments. There is another argument 
why the information paradox shouldn’t be related to quantum gravity at all. In US a bound 
for shear viscosity over entropy density was proposed, r\f s > 1 /( 47 r), and in |2T it was shown 
that the bound is saturated by black holes in gravity duals. But in RHIC collisions, 77 /s is 
close to the limit value [221; a fact that we used in |3J to support our assertion that RHIC 
fireballs are gravity dual black holes living on the IR brane. The gravity dual saturation of 
the bound was derived using apparently quantum gravity arguments for the thermodynamics 
of black holes, yet the bound itself is independent of the Planck scale Mpi (this being the 
reason why it is possible to have the same bound saturated at RHIC). That would suggest 
that the thermodynamics of black holes might not have anything to do specifically with 
quantum gravity, but might just come out of usual held theory. 

In section 7, we observed that the correct effective description in the real (QCD) vacuum 
is given by the “dual superconductor” picture of’t Hooft. The dual of the maximal abelian 
subgroup 17(1) x 17(1) of 517(3) is in a type II superconductor phase, the dual U(l) gauge 
fields being massive due to a Higgs held, except for hux tubes between dual monopole- 
antimonopole pairs representing mesons, with baryons being some spherically symmetric 
version of this picture. Thus the dual effective picture a la Seiberg-Witten for the nucleons 
in the usual vacuum should involve the U(l) gauge helds, the Higgs making them massive 
and the scalar SU(2) pion helds. Assuming that we can neglect the dynamics of the Higgs 
and postulating the DBI high energy form of the action based on the gravity dual as well 
as Heisenberg’s model, we wrote down an effective action for QCD in the usual vacuum, as 
given in O' It is a DBI action for the U(l) massive gauge held coupled to SU(2) light 
pions. The action has the nonsingular Blon solution, however, as soon as we add a scalar 
perturbation at infinity, the solution will be singular at a hnite radius, having a “horizon”, 
with hnite scalar F, but inhnite F’. As a result, there doesn’t seem to be any topological 
solution to the action, only solutions with horizons. However, as argued, the DBI action 
must be modified where the energy density of the U(l) gauge held diverges, and thus we 
argued that there will be a solution to the modified action, that we dubbed skyrBIon. 

We argued that the topological charge B is the baryon charge, as usual, with B=+l being 
baryons, B=-l antibaryons and higher B solutions, if present, representing states unstable 
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against decay into B=1 states because of energy considerations. The objects of charge +1 
and -1 under the massive U(l) should represent the n and the p (the n and the p have 
different quark content, which will therefore have different gluon interactions). We argued 
that the potential between two nucleons should have a minimum at an rc of the order 
°f 1/A qcd 0.2/m. This feature is enough to guarantee that a nucleus can form from 
individual nucleons, in a liquid drop-like model. Given that Blons form NaCl-like crystals 
of alternating charges, we argued that the nuclear energy for skyrBIon nuclei should be 
minimized for N=Z, giving a qualitative agreement with the Bethe-Weizsacker formula. 

For the high energy (high s, fixed t) scattering of skyrBIons we could find qualitative 
agreement with the gravity dual picture, by identifying |VT(r)| with £{r) = y/sa(r), but 
the particular small power laws obtained in the gravity dual cannot be deduced just from 
the effective DBI action. In the high energy scattering of two baryon-like solutions of the 
U(l) D-brane toy model one will create a catenoid= scalar fireball of temperature that 
goes like T oc M -1 / 3 , where M is the total energy of the collision. In this regime, the 
fireball does not behave yet like a black hole, except for having a temperature. In the 
case of skyrBIon collisions, the created fireball solution will be similar. We have argued 
for the existence of a “soft Pomeron scale” and a Froissart onset scale, as we derived from 
the gravity dual picture, but we can’t calculate them, except to say that both should be 
larger than Mj>/ A.q CD = NAqcd■ The Froissart regime is associated with production of a 
scalar fireball of size much larger than 1/m^, and the temperature of the fireball should be 
proportional to m n in this asymptotic regime, as argued in [Hj from a gravity dual point of 
view. Then the produced scalar fireball is directly mapped to a gravity dual black hole living 
on the 4 dimensional IR brane. 

Thus as we advocated in [S], we have shown that we can describe the gravity dual picture 
for high energy scattering completely in terms of field theory, but we have seen the limitations 
of field theory in terms of calculability. The thermal property of black holes in the gravity 
dual is easily understood as the thermal property of “horizons” in the effective pion field. 
The colliding pion field shockwaves are seen as being just boosted versions of nucleons. The 
nucleons are described as solutions of an effective action, and lead to good qualitative features 
for the description of nuclei. 
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